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A unified normal mode approach to dynamic tides and its 
application to rotating Sun-like stars 
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ABSTRACT 

We determine the response of a uniformly rotating star to tidal perturbations due to a com- 
panion. General periodic orbits and parabolic flybys are considered. We evaluate energy and 
angular momentum exchange rates as a sum of contributions from normal modes allowing for 
dissipative processes. We consider the case when the response is dominated by the contribu- 
tion of an identifiable regular spectrum of low frequency modes, such as rotationally modified 
gravity modes. We evaluate this response in the limit of very weak dissipation, where indi- 
vidual resonances can be significant and also when dissipative effects are strong enough to 
prevent wave reflection from the neighbourhood of either the stellar surface or stellar centre, 
making radiation conditions more appropriate. The former situation may apply to Sun-like 
stars with radiative cores and convective envelopes and the latter to more massive stars with 
convective cores and radiative envelopes. We provide general expressions for transfer of en- 
ergy and angular momentum that can be applied to an orbit with any eccentricity. 

Detailed calculations require knowledge of the mode spectrum and evaluation of the 
mode overlap integrals that measure the strength of the tidal interaction. These are evalu- 
ated for Sun-like stars in the slow rotation regime where centrifugal distortion is neglected in 
the equilibrium and the traditional approximation is made for the normal modes. We use both 
a WKBJ procedure and direct numerical evaluation which are found to be in good agreement 
for regimes of interest. The former is used to provide expressions for the mode spectrum and 
overlap integrals as a function of mode frequency and stellar rotation rate. These can be used 
to find the tidal energy and angular momentum exchange rates and hence the orbital evolution. 

Finally we use our formalism to determine the evolution time scales for an object, in an 
orbit of small eccentricity, around a Sun-like star in which the tidal response is assumed to oc- 
cur Systems with either no rotation or synchronous rotation are considered. Only rotationally 
modified gravity modes are taken into account under the assumption that wave dissipation oc- 
curs close to the stellar centre. It is noted that inertial waves excited in the convective envelope 
may produce a comparable amount of tidal dissipation in the latter case for sufficiently large 
orbital periods. 

Key words: hydrodynamics - celestial mechanics - planetary systems: formation, planet -star 
interactions, stars; binaries: close, rotation, oscillations, solar-type 



1 INTRODUCTION 

The tidal interaction between two bodies is an important problem 
in astrophysics. It plays a role in close binary stars leading to syn- 
chronisation and orbital circularisation (eg. Zahn 1977, Hut 1981), 
as well as stars undergoing a close encounter that may lead to a tidal 
capture (eg. Press & Teukolsky 1977). In recent years the discovery 
of many extrasolar planets orbiting in close proximity to their cen- 
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tral stars has highlighted another situation where tidal interactions 
may play an important role in the formation and evolution of the 
systems (eg. Barker & Ogilvie 2009). In particular hot Jupiters may 
be formed by a tidal capture into a highly eccentric orbit followed 
by orbital circularisation (Weidenschilling & Marzari 1996, Rasio 
& Ford 1996, Nagasawa, Ida & Bessho 2008). During this process 
tidal dissipation both in the planet and the star can be significant 
(Ivanov & Papaloizou 2004, 2008, 2011). 

The general problem of determining the tidal evolution of a 
system involves finding the response of the rotating stars and/or 
planets due to the gravitational perturbation of a companion. This 
results in the excitation of normal modes and the exchange of en- 
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ergy and angular momentum with the orbit leading to its evolu- 
tion. This may be technically challenging on account of the nor- 
mal mode spectrum being dense and/or continuous (Papaloizou & 
Pringle 1982). Critical latitude phenomena and wave attractors may 
also be involved in some cases (Ogilvie & Lin 2004). Multidimen- 
sional tidal response calculations need to be undertaken for bodies 
in a variety of evolutionary states and rates of rotation. Accordingly 
it is useful to formulate methods for calculating the tidal interac- 
tion that involve the evaluation of manageable analytic expressions 
involving functions of at most one variable that can be readily tab- 
ulated. 

In this paper we develop general procedures for expressing 
tidal energy and angular momentum exchange rates in simplified 
terms that are applicable to the tidal response that arises from an 
identifiable regular spectmm of low frequency modes, comprised 
of for example rotationally modified gravity modes, for rotating 
stars with realistic structure. These are likely to give the dominant 
tidal response in bodies with stratification, where the tidal forcing 
frequencies significantly exceed the inverse of the convective time 
scale associated with any convection zone, making the associated 
turbulent viscosity inefficient. The waves associated with the rota- 
tionally modified gravity modes may either produce a long lasting 
resonant response through the so-called 'mode locking' mechanism 
(eg. Witte & Savonije 1999, 2002) , or be dissipated, either near the 
centre in Sun-like stars with a radiative core and convective enve- 
lope (eg. Barker 201 1), or near the surface in massive stars with a 
convective core and radiative envelope (eg. Zahn 1977). 

When decomposition of free stellar pulsations over normal 
modes is possible, and the damping mechanism of normal modes is 
specified, we give general expressions from which the energy and 
angular momentum exchange may be calculated, in principal, for 
both bound orbits of arbitrary eccentricity and parabolic encounters 
that could lead to tidal captures with an arbitrary inclination of stel- 
lar rotation axis with respect to orbital plane. These are expressed 
in terms of contributions from individual normal modes. In general, 
they require evaluation of coefficients characterising the represen- 
tation of the tidal potential as either a Fourier series or a Fourier 
transform, the spectrum of normal modes, normal mode damping 
rates and so called overlap integrals that measure the strength of 
interaction. It is important to note that in certain regimes of tidal in- 
teraction, such as energy and angular momentum exchange through 
parabolic encounters, or orbital evolution in the regime of so-called 
moderately strong viscosity (eg. Zahn 1977, Goodman & Dickson 
1998), the resulting expressions do not depend on the details of the 
mode damping, which is poorly known in many situations, see also 
the discussion of these regimes below. 

The overlap integrals determine the coupling of mode eigen- 
functions with the tidal potential, they thus play a central role in our 
approach to tidal problem. Note that although we use the definition 
of overlap integral first introduced by Press & Teukolsky (1977) 
for non-rotating stars and later generalised by Ivanov & Papaloizou 
(2005, 2007) for rigidly rotating stars, they are directly related to 
the so-called tidal resonance coefficients first introduced by Cowl- 
ing (1941) and studied later in a number of papers (eg. Zahn 1970, 
Rocca 1987, Smeyers, Willems & Van Hoolst 1998). In order to 
calculate the overlap integrals one should be able to evaluate inte- 
grals with highly oscillatory integrands over the bulk of the star. 
This is rather difficult to achieve both numerically and analytically. 
In such a situation, numerical and analytical methods should sup- 
plement each other. In this paper we give an analytic method for 
their evaluation, together with a calculation of the relevant eigen- 
spectrum and normal modes, valid for Sun-like stars having an in- 



ner radiative region and a convective envelope, that is based on a 
WKBJ scheme. In an accompanying paper we also provide values 
of overlap integrals for a range of eigenfrequencies in a number of 
models of more massive stars of different masses and ages. 

It is quite important to stress that the calculation of the rel- 
evant normal mode spectra for rotating stars is complex, because 
unlike the situation for a spherical body, separation of variables is 
not in general possible. However, this is rectified if the well known 
traditional approximation is used as is done below. In the tradi- 
tional approximation only the radial component of the angular ve- 
locity is retained (Unno et al. 1989). This expected to be a good 
approximation for rotationally modified gravity modes for which 
motions are mainly horizontal. It has accordingly been adopted by 
many authors. For example, it was used in an asymptotic treatment 
of the normal mode problem by Berthomieu et al. (1978) and in di- 
rect tidal response calculations (Papaloizou & Savonije 1997, Witte 
& Savonije 1999, 2002 and Ogilvie & Lin 2004). However, it is 
worth mentioning that the existence of regular spectra of rotation- 
ally modified gravity modes is not dependent on the use of the tra- 
ditional approximation. Such spectra are expected to be associated 
with the radiative cores of Sun-like stars and the radiative envelopes 
of more massive stars. A spectrum of rotationally modified gravity 
modes in the case of a massive slowly rotating star was identified 
by Savonije et al. (1995). In this case the traditional approximation 
was found to be a good one as long as their eigenfrequencies placed 
them outside the inertial regime ( the absolute value of the eigen- 
frequency exceeded twice the stellar angular velocity). Similarly 
Dintrans & Rieutord (2000) found a rotationally modified gravity 
mode spectrum outside the inertial regime for a l.SA/^ star with a 
convective envelope without making use of the traditional approx- 
imation. Note too that the order of the I = 2 g mode with eigen- 
frequency corresponding to the critical break up rotation frequency 
for a Sun-like star is typically ^ 10. Accordingly a low frequency 
asymptotic analysis should be useful for stars in a state of modest 
rotation, at less than about one fifth of the critical rate, and excited 
modes with comparable eigenfrequencies. 

Finally, we demonstrate how our approach can be applied to 
the calculation of the orbital circularisation and decay time scales 
for objects orbiting Sun-like stars under the assumption of wave 
dissipation near the centre and small orbital eccentricity. 

The plan of this study is similar to that of our previous study 
of inertial waves in rotating bodies undertaken in Ivanov & Pa- 
paloizou (201 1) and Papaloizou & Ivanov (201 1). In this paper we 
mainly concentrate on an analytic approach to the problem while in 
an accompanying paper we consider the application of numerical 
methods, extending our calculations to evaluate the tidal response 
of stellar models appropriate to stars more massive than the sun. 

In section [2] we give the basic equations governing the tidal 
response of a rotating star to tidal forcing. We then use them to cal- 
culate the linear response to harmonic forcing for a particular az- 
imuthal mode number in section|3] To do this we apply an extension 
of the method of Ivanov & Papaloizou (2007), hereafter IP7, which 
provides a solution in terms of a decomposition in terms of normal 
modes, that allows us to incorporate dissipative forces. Specifically 
we consider viscous forces for which the viscosity coefficient may 
vary with location but the approach may be generalised to deal with 
more general forms of dissipation. 

We focus on the case where there is a regular spectrum of nor- 
mal modes that can be dense, such as a low frequency gravity mode 
spectrum and consider its response. In particular we calculate the 
rate of transfer of energy and angular momentum between the star 
and the orbit, expressing the results in terms of overlap integrals 
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that measure the contribution associated with each normal mode in 
teiTns of the projection of a related eigenfunction onto the tidal po- 
tential. We then customise the expressions for the tidally induced 
energy and angular momentum exchange rates to apply to the case 
where the response is determined by a dense regular set of normal 
modes, such that the mode decay rates are small compared to their 
natural frequencies, but not small compared to the separation be- 
tween adjacent mode frequencies. 

As we demonstrate in an appendix for slow rotation under the 
traditional approximation, this case is expected to correspond to 
the situation where waves associated with the modes propagate to 
localised boundary regions, which may be either close to the stel- 
lar centre or surface, where they are dissipated, and we call it the 
moderately viscous case. This is in contrast to the situation where 
dissipation is very weak, and the interaction is dominated by strong 
normal mode resonances that are effective only in very narrow fre- 
quency intervals (eg. Savonije et al. 1995, Witte & Savonije 1999, 
2002). The expressions we obtain can be used to determine the or- 
bital evolution of a binary system, in a general periodic orbit which 
may be near circular, or undergoing encounters in a near parabolic 
orbit, induced by tides. Knowledge of the adiabatic normal mode 
spectrum and the overlap integrals corresponding to each mode, in 
the low frequency regime, is required for their evaluation in the 
moderately viscous case, the results then being independent of the 
details of the dissipation process. 

In section |4] we go on to discuss the linear adiabatic normal 
mode problem and the low frequency spectrum provided by rota- 
tionally modified gravity modes. We consider the low frequency 
asymptotic regime under the traditional approximation in which 
only the radial component of the angular velocity is retained. 
When this approximation is made, separation of variables is pos- 
sible. Its efficacy in the stellar case was first implicitly realised by 
Berthomieu et al. (1978), through the application of an asymptotic 
procedure, that effectively led to it. We develop expressions for the 
overlap integrals in the low frequency limit that can be evaluated 
as one dimensional integrals involving the mode angular and radial 
eigenfunctions. 

In section|5]we perform direct numerical calculations of mode 
eigenfrequencies, eigenfunctions and their associated overlap in- 
tegrals. Rotationally modified gravity modes under the traditional 
approximation are found for models of solar-type stars. Models 
representing both the present day sun and a young sun with age 
1.66 X lO^yr. are considered. We also give power fits to the de- 
pendence of the square of the Brunt - Vaisala frequency inside the 
radiative core, on the distance to the boundary between the interior 
radiative zone and the convective envelope (see also Barker 201 1). 
This dependence was found to have a significant effect on proper- 
ties of the normal modes. 

In section |6] we calculate eigenfrequencies, eigenfunctions 
and overlap integrals for high order adiabatic rotationally modified 
gravity modes for models of solar-type stars under the traditional 
approximation following a WKB J approach that incorporates some 
corrections to the formalism of Berthomieu et al. (1978). We also 
allow for an arbitrary power law dependence of the Brunt - Vaisalla 
frequency in the radiative interior on the distance to the convective 
envelope boundary. We provide expressions for both eigenfunc- 
tions and eigenvalues which are found to provide a good approxi- 
mation to the corresponding numerically determined quantities for 
the stellar models considered. 

We use these solutions to evaluate the overlap integrals, pro- 
viding analytic expressions for them, which together with some 
quantities represented graphically, readily enable their evaluation 



as a function of mode frequency and rotation rate of the star. These 
are in very good agreement with the numerically determined quan- 
tities for mode frequencies less than the critical rotation frequency. 

Knowledge of the overlap integrals and mode spectrum en- 
ables calculation of tidal energy and angular momentum exchange 
rates with a general orbit. In section |7] we apply our formalism to 
the calculation of the tidal evolution timescales for a binary with 
a small eccentricity and with aligned orbital and spin angular mo- 
menta. One of the components is assumed to be a Sun-like star 
in which modified gravity modes are excited and propagate to the 
centre where they are assumed to be dissipated corresponding to the 
moderately large viscosity regime. The other component is treated 
as a point mass and a large range of mass ratios is considered. We 
give estimates for the circularisation time scale when the star is 
non rotating and when it rotates synchronously with the orbit. In 
the former case we give the time scale for the semi-major axis to 
decrease. 

Finally in section[8]we summarize and discuss our results. 



2 BASIC EQUATIONS 

2.1 Coordinate system and notation 

The basic definitions and notation adopted in this paper closely fol- 
low IP7. We use either cylindrical coordinate system {zu, 4>, z) or 
spherical (r, </>, 9) coordinate systems with origin at the centre of 
mass of the star. When viewed in an inertial frame with this origin, 
the unperturbed star is assumed to rotate uniformly about the z axis 
with angular velocity Q. Accordingly we adopt the rotating frame 
in which the unperturbed star appears at rest as our reference frame. 

In this paper the fundamental quantity describing perturba- 
tions of a rotating star will be the Lagrangian displacement vector 
^. By definition, this is a real quantity. We represent the displace- 
ment vector as well as any other perturbation quantity in terms of 
a Fourier series in the azimuthal angle </!>, and either in terms of a 
Fourier series in the time, t, or a Fourier integral with respect to 
t, depending respectively on whether we consider tidal interactions 
in a binary having a periodic orbit, or the orbit is assumed to be 
parabolic. 

When a discrete Fourier series in time is appropriate, for a real 
quantity, Q, we write 



m,k 



{Qm,k exp{—iujm,kt + im4>) + cc) , 



(1) 



where cc denotes the complex conjugate, ojm.fc ~ kflorb — mQ, 
Qorb is the orbital frequency, Q is the angular velocity of the star 
and m and k are integers, with only positive values of m included 
in the summation as in IP7. Note when the orbital plane is perpen- 
dicular to the rotational axis of the star, and only the quadrupole 
component of the tidal potential, that is dominant at large separa- 
tions is included, terms with 7n = 1 are absent in the summation. It 
then suffices to consider jn = 0, 2. In addition, the reality of Q im- 
plies that we should identify Q _„,-/.- = Qm,k- When forming the 
sum ([TJ, we ensure that terms are not repeated. For example, when 
m = Q, terms for which uirn.,k changes sign are identical, they are 
accordingly combined. After this is done, forcing for only one sign 
ofujrn,k has to be considered . 

For the case involving a Fourier transform, we write 



Q =y I exp(im 



doQm exp(— iat) + cc 



(2) 
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where the Fourier transform is indicated by a tilde. This require- 
ment that Q is realleads to the identification (5_m(—o-) — (5m(o")- 
The inner products of two complex scalars Yi, Y2, and two complex 
vectors r]^ and 773 are respectively defined through 



and (3) 

(4) 



where (771* -772) is the scalar product, p is the density, and * stands 
for the complex conjugate. 



2.2 Equations of motion 

In the rotating frame the linearised Navier-Stokes equations take 
the form (see IP7) 



(YiiVa) = / zudujdzp{Y^Y2) 
(»7il»72) = wdujdzp(ril-ri2), 



^ + 2nx^ + C^ = -VvI/-/,, 



(5) 



where the dot indicates partial differentiation with respect to time 
and C is an integro-differential operator describing the action of 
gravity and pressure forces. For conservative boundaiy conditions, 
it is self-adjoint when the normal product (|4j is adopted, and it can 
be considered to be a positive operator under certain conditions, 
which are assumed to be fulfilled in this paper. Its general explicit 
form is not important for our purposes and can be found in e.g. 
Lynden-Bell & Ostriker 1967. 

Retaining only the assumed dominant quadrupole component, 
the tidal potential, ^ can be represented in the form 



* = r 



rn,,k 



.,fce 



'KT + cc) 



(6) 



where Y^J^ are the spherical functions. The coefficients A^.k are 
given in Appendix A for the case of coplanar orbit with small ec- 
centricity. For orbits having an arbitrary value for the eccentricity, 
these coefficients can be expressed in terms of so-called Hansen co- 
efficients discussed eg. in Witte & Savonije (1999). However, note 
that the representation of the tidal potential in terms of a Fourier 
series ([6f is not convenient for a highly eccentric orbit. In that case 
it is more appropriate to adopt the tidal potential for a parabolic or- 
bit, and its Fourier transform. The corresponding expressions can 
be found in e.g. Ivanov & Papaloizou (201 1). These authors as well 
as Lai & Wu 2006 also give expressions, which can be used to rep- 
resent the tidal potential in a coordinate system that is rotated with 
respect to the one for which the z axis is normal to the orbital plane. 
These expressions can be used when the rotation axis of the star is 
not perpendicular to the orbital plane. 

The viscous force /„ = cr^, where the operator a is defined 
through its action on a vector rj as 



(Jafi = Va,l3 + V0,a 



-r]-t,-tSc0 (7) 



p 

where v is kinematic viscosity, a comma stands for differentiation 
over Cartesian coordinates and summation over repeating indices 
is implied. Note that we reserve the indices a, /3 and 7 to indi- 
cate components of vectors and tensors in order to distinguish them 
from indices used in the summation of series. When these indices 
are repeated, summation is implied from now on (the summation 
convention). 

We remark that we have considered a standard form of vis- 
cous dissipation with the possibility that the viscosity coefficient 
depends on position. This leads to some simplifying features in the 



discussion of dissipation below on account of the viscous stress ten- 
sor being symmetric. However, a corresponding analysis can also 
be undertaken for more complicated forms of dissipation, including 
radiative damping, by introducing relevant adjoint operators. How- 
ever, as the final results we use do not depend on the details of the 
dissipation beyond involving decay rates associated with individual 
noiTnal modes we shall work only with the standard formulation for 
viscous dissipation described above. 



3 THE RESPONSE TO TIDAL FORCING DUE TO 
COMPANIONS IN PERIODIC AND PARABOLIC 
ORBITS 

In this section we calculate the tidal response to tidal forcing asso- 
ciated with a general periodic orbit. Later we extend the discussion 
to the limiting case where the orbit becomes parabolic (see IP7). 
We begin by writing the forcing potential (|6) as a Fourier series of 
the form ([T}. The Fourier coefficients, iprn,k, are such that 



(V*)™,fc 



A 



m,k 



2ti 



-imcpyj / 2-tr'rn 



^{rYD, 



(8) 



where Y^^ indicates the usual spherical harmonic with indices m 
and I. 

As was shown by Papaloizou & Ivanov 2005 and IP7, after 
taking the Fourier transform of equation (|5], one can obtain equa- 
tions for determining the Fourier coefficients of the Lagrangian dis- 
placement ^^ J. . These can be organised into a convenient form as 
described below. Once determined for periodic orbits, the ^^ j. can 
then be used to calculate the orbital energy and angular momentum 
exchange rates. 



3.1 Determination of the Fourier components of the 
Lagrangian displacement 

We begin by introducing a six-dimensional vector Zm,k 
(Zi,Z2)^ where 



(9) 



Zi=ijJm,k$,^^k^ and Z2 = C^/^^ 

Equation ([Sjl then follows from the pair of equations derived from 

^m,kZm,k = 'HZm,k — i^m,kT^ Zm,k + S, (10) 

where the operator "H has a matrix structure 



■H 



ei^ 



(11) 



with the operator C^ " being defined through the condition: C = 
C^' C^' , with any ambiguity removed by requirement that it is a 
positive operator ( see eg. Wouk 1966). In addition the operator 



Bi = -2iQ, X i. 

The operator "D is defined through 

T^Zm,k = {fvx,k^iv-i,k) > 

where 



/--I fe = f^-^i and j^^ fc = 0. 



(12) 



(13) 



(14) 



As we have performed a Fourier decomposition in terms of 
m, the action of the operator O being any one of B,C or ct on the 
Fourier component ^^ ). is obtained through 
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C^m.fc = exp(-i?n(?!))0(^„_s^ exp(im(;/))). 



(15) 



It readily follows that B is self-adjoint as are C and C^ " (IP7) . 
From the definition (|7} it also follows that a is self-adjoint. The six 
dimensional source vector S = {Si, S2)''" where 



Si = (V*)™,fe and 52 



0. 



(16) 



Equation JlOb can be used to define an eigenvalue problem in 
the form 



X^jZ^ = (■H-iujm.k1>)Z' 



(17) 



Here Xj is the generally complex eigenvalue and the corresponding 
eigenfunction is Z-' = {Z\j, Z2j)'^ , the components being related 
to the corresponding displacement by 

and Z2i=C^''^'- 



Zi, =A,^^., and Z2,=C''^£.,. 
Adopting the inner product defined through 

< X\Z >^ zudvjdzp{X*i,X2)iZi,Z2 



(18) 



(19) 



it follows from equations dllt and ( I13t that the operators "H and "D 
are self-adjoint. Hence it follows that the operator 'H + iu)m,k'D is 
the adjoint operator associated with the eigenvalue problem l ll7b . 
The eigenvalues will accordingly be eigenvalues X*j with corre- 
sponding eigenfunctions X-' . 

In addition the eigenfunctions Z-' and X-' are biorthogonal 
with respect to the adopted inner product so that 



< P\Z' >= X*MriMd + (r/,|CC,) = S,,N, 



(20) 



where rj- denotes the Lagrangian displacement associated with the 
eigenfunctions of the adjoint problem defined through ( 118) with Z 
being replaced by X. The norm is 



iV, =<X^1^^- >=|A,|2(r,,|e,) + (r7,|C^,)• 



(2I) 



We remark that although we have not assumed the viscosity is small 
up to now, when it is neglected ^ and r; ■ coincide. 

3.2 Solution for the forced response 

We look for solution of dlOt as a decomposition over the eigenvec- 
tors Z' : 



-'vn-.k 



Y.^^^' 



(22) 



Substituting l l22t in MQ\ , taking the inner product of the result with 
the eigenvector X' and using Ml\ and ( |20t we get 



ai{uim,k - Xi)Ni =<X'|5'> 



(23) 



To proceed further we remark that an expression for the damp- 
ing rate of the mode with eigenvalue Xj can be found by taking the 
inner product of ( II7t with Z-' . Using the fact that T-i is self-adjoint, 
we find that 



Tm{Xj) < Z^\Z^ >= -oj^,k < Z^-DZ^ >, 



(24) 



where Xm denotes the imaginary part. Setting Im{Xj) — —co^^kj 
and making use of JI8t . ( I13t and Q we obtain 






(25) 



where the superscript j indicates that (7^„ is evaluated assuming 
that r) as used in equation l|7} is replaced by ^ and 



Using the above we write A, — 



Ov,kj 



(26) 
where oj, is the real 



part of Xj. We remark that for convenience in the above discus- 
sion we have suppressed the index k for Xj , ujj and ^ . Then from 
equation l l23t we find 

_ Xjbj^k 



Nj(bJm,k - t^j + il^v.kj) ' 

where 

Sj.k = (T/j|V*™,fc), 

Now we use equation ( 122b to obtain 



>m,fc 



■^jSj.k 



^-^ ^m,kNj{uJm.,k — ^J + i^v,kj) 
3 



(27) 



(28) 



(29) 



3.3 Energy transfer rate due to tidal forcing 

The canonical energy associated with the perturbations is given by 
a standard expression 



E, 



d'^^Pl^(\k.? +$,-C$,) which leads to 



(S,) = 27r^«,(^™,fel^^,fc) + (^™,fe|C^,„.J}, (30) 

Tn , k 

where we integrate over the volume of the star Recall that ^ is 
real. The last expression, being a time average, displays explicitly 
contributions corresponding to different azimuthal numbers m and 
forcing frequencies ujm.k- The contribution of an equal response 
to the complex conjugate of the forcing potential in (|6]l has also 
been included. The rate of change of Ec with time as a result of the 
action of dissipative forces follows from (|5) and ( |30t as 



Ec ^ - I <fxp{^ ■ f^) = -- / (fxpu[aaB<7a,B] 



(31) 



where we use (|7} with r] replaced by ^ and integrate by parts. We 
remark that when, as here a sum of steady periodic responses is 
excited by a forcing potential, the rate of energy loss resulting from 
dissipative forces is balanced by input from the forcing potential. 
We now go on to obtain an expression for Ec in the form of a sum 
of contributions corresponding to different m and uJm.k- 

To do this we substitute ( |29l l in the appropriate form of ([T) and 
the result into ( I31t . An average of the resulting expression over a 
time period that is large in comparison to the periods of eigenmodes 
that contribute significantly to the sum is then taken. We thus obtain 



Ec 



where 






J j' '^j'^^ ^' 



3 ,k 



NjN*,D. - -, 



\ a 



I3l'^a0. 



+ CC 



D 



k,3,3 



= {^m,k — ^j + i^v,kj){^m,k — ^,' — *^^.fc,' ) 



(32) 



(33) 



When the mode energy does not grow with time on average, 
and, accordingly, the stationarity condition is implied as in the bulk 
of the papeo the rate of energy dissipation by the action of viscosity 
must be negative of the work done by tidal forces, Ec ~ —Et, 
where Et = — 27r(^jV^), see equation ([S} and averaging over 



^ We briefly consider the situation when the stationarity condition may be 
invalid in the end of Discussion. 
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a time larger than a characteristic time scales associated with the 
forcing frequencies is implied. From equations |[T), (|9) and M2\ it 
is easy to see that 



Et = 27ri ^ ^m,fc(Qj((V*)„,fcl^,„ ,.) - cc) 



(34) 



7n,k,j 



Note that when viscosity is small the complex conjugate of 
(^m k\V^'^)m,i<:) IS approximately given by equation ( l28t . 



3.3. 1 The case of very small viscosity and near normal mode 
resonances 

The expression ( |32t can be further simplified by assuming that a 
strong resonance exists in the vicinity of a particular forcing fre- 
quency, uJm,k, on account of very weak dissipation. If we assume 
as a consequence of this that, for a given forcing frequency, the 
dominant contribution to the sum comes from a particular mode 
that is nearly resonant, only the diagonal terms in the sum over j 
and j need to be retained as they provide this contribution. But 
note that the appropriate damping rates due to viscosity should be 
sufficiently small, as indicated further below. Then we have 



E^ 



m.k,j 






3,k\ 



m^D 



k.j.j 






'afi\'^afi. 



\SjA' 



m.k,j 

where we have used equation ( l25t to obtain the last equaUty. 



(35) 



3.3.2 The case of moderately large viscosity 

In order to make the above procedure valid, the effect of viscosity 
should be small. But note that we shall be interested in the case 
when the viscous decay rate of a mode, although being small com- 
pared to its frequency, may be large compared to the frequency sep- 
aration between modes. This situation is less clear as many modes 
could interact and treatment of viscosity by straightforward pertur- 
bation theory is not possible. However, taking the inner product of 
equation l ll7t with Z^ and subtracting the complex conjugate of 
the corresponding equation with j and j' interchanged, we obtain 



(A, - \*,) < Z' \Z' >= -2ioj„^.k < Z' \VZ' > 
which leads to the following expression 



(36) 



i(A.A*'(^,.|^,) + (^,.|C^,)) 



(37) 



(A,- A*,) 

From this we can argue that, when the mode separation fre- 
quencies are less than the characteristic damping rate (assumed to 
be comparable for neighbouring modes), the contribution of off di- 
agonal terms, as compared to diagonal terms in ( 132) is measured 
by the magnitude of an overlap expression between the displace- 
ments of the two modes that approaches zero in the limit of zero 
viscosity when j ^ j' . Provided such overlaps remain relatively 
small, it should be possible to neglect the off diagonal terms in the 
sum. This would be satisfied if the form of the eigenfunctions is 
preserved when the level of viscosity of interest is introduced. This 
is physically reasonable when the damping rate is small compared 
to the actual mode frequency and viscosity is important only in 
localised regions of the star/planet where the mode wavelength is 
small, that for example, could be modelled using a boundary layer 



treatment, in such a way that the form of the normal modes of in- 
terest is preserved in most of the star/planet. Such a situation may 
occur for stars with radiative envelopes when g modes are effec- 
tively damped only near the surface. We shall assume it to be the 
case from now on so that ( I35t applies and inviscid modes can be 
used to evaluate the terms in the summation. 



3.4 Angular momentum transfer rate and energy transfer 
rate as viewed in the inertial frame 

Equation ( 135b gives the transfer rate of canonical energy defined in 
the rotating frame. The angular momentum transfer rate, Lc, and 
the energy transfer rate as defined in the inertial frame, Ei, follow 
by using the fact that for forcing with a harmonically varying po- 
tential with pattern speed ujm.k/'m, the rates of change of angular 
momentum and canonical energy are related by Lc — mEc/u)m,k- 
In addition, for a harmonically varying forcing potential, the rate 
of change of energy in the inertial frame is related to the rates of 
change of canonical energy and angular momentum through 



El — Ec + flLc — Ec{ljJ,n± + mQ,)/ujm,k. 



(38) 



Applying the above results separately to the contributions to the 
sum in equation ( I35t arising from distinct forcing frequencies, we 
obtain 



-47r 



E 

m,k,j 



Ei = -47r V" 

m,,k,j 



muj^,kj\Xj\^J^j\Sj,k\'^ 

^m.,k\Nj\^Dk,j,j 

^»^,fcj|AjP(^m,fc + mQ)Afj\Sj,k\ 

LJm,k\Nj\^Dk,j,j 



(39) 



(40) 



Note that in the limiting case of a small stellar angular velocity di- 
rected peipendicular to the orbital plane, and under the assumption 
that the effect of rotation can be dealt with by simply Doppler shift- 
ing forcing frequencies, analogous equations have been derived by, 
e.g. Kumar etal. (1995). 



3.5 A simplification in the case of dense spectrum of normal 
modes 

The expressions ( I35t ( 139) and ( 140) can be significantly simplified 
in the case when the spectrum of eigenmodes contributing most 
significantly to the sums becomes dense with a typical distance be- 
tween two neighbouring modes \ujj+i — ijJj\ approaching a smooth 
function that we can identify with \dhjj/dj\ <C \ujj\ for j — > oo. 
An example of this situation occurs for non rotating cool stars with 
radiative regions that support internal gravity modes. These have 
the required property in the high order/low frequency limit. In that 
case \dujj/dj\ — >■ \ijJj/j\ for j — >■ cxa. In general, a high degree of 
regularity of the spectrum is required for it to have this property. 

As above we suppose that the main contribution to the series 
comes from modes having near resonance with the external forcing 



frequencies. These are such that ljj 



fc making the magnitude 



of the denominator, Dkjj , a small quantity for such modes. Let us 
assume that for a particular forcing frequency, (jJm,k, the resonance 
condition is most closely satisfied for a particular mode with index 
j — jo, (which we recall will vary with m and k). Accordingly we 
have 

ujj„ = aj„,fe + Atjjo, (41) 

where the offset, Acj^g, is such that lAajj,,] < \dujjg/djo\. We 
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approximate values of neighbouring mode frequencies by repre- 
senting ujj as a Taylor series in j — jo up to the linear teiTn, thus 
Ljj ~ cjjo + {dbjjg/djo)l, where ^ = j — joQ This is then substi- 
tuted in Dk,j,j and the result used in equations ( I35t ( 139) and ( I40t . 
In addition quantities other than Dk,j,j in eg. | |35I ) are assumed to 
depend smoothly on j. As the dominant contribution to the sums 
comes from values of j close to jo , such quantities can be evalu- 
ated for j = jo and Uj = ujm,k and taken out of sums over j. 
In addition as the viscous decay rates are assumed small we set 
Xj = LJm,k- Following the discussion in section [3.3.21 we assume 
the normal mode spectrum is approximately independent of dissi- 
pation so that the inviscid case may be adopted, we set Nj and A/'j 
equal to A/'jo = N jo- 
in this way from ( I35t we obtain 



E, = 



Ak 



"m.k 



m.k 



'jo.fcl 



N„ 



where 



^ {i + sy 



and we formally extend the sum over the index [l) to ±00, 



5 = 



Ao;,, 



dujja/djo 



and K = 



"■^.fejo 



dujo/djo 



(42) 



(43) 



The factor Ak can be evaluated by standard complex variable tech- 
niques to give 

■K sinh TTK cosh ttk 



A. 



K(sinh^ TTK + sin"^ tvS) 



(44) 



3.5.1 The case of moderately large viscosity or an averaged 
dissipation 

Let us consider the case of large k ^ 1 which occurs when 
the characteristic viscous decay rate of a resonant mode is much 
larger than the frequency difference between the relevant mode and 
a neighbouring mode. On the other hand the viscous decay time 
should also be much longer than the inverse mode frequency so 
that strong resonances with one particular mode are still possible. 

As discussed by Goodman & Dickson (1998), physically this 
situation occurs when the decay time of a mode due to viscosity is 
much smaller than the time for wave propagation with the group 
velocity, associated with the potentially resonant modes, which is 
identified as |dajjo/d/o|~^- To see that the latter identification is 
reasonable, suppose that mode j is associated with an appropriately 
defined mean wavenumber kj which, for large j, asymptotically 
approaches j/A, where A is a length scale that is characteristically 
of the order of the size of the region containing the modes. We have 
\duiJQ /djo I = \du)j /dkj j j=jq /A. It is seen that the right hand side 
of the above is the inverse of the time to propagate across A with 
the notional group velocity \dujj/dkj |j=jq . 

In the limit k 3> 1, Ak ~ Tr/fi: and the expression for the 
energy transfer 



Ec = -47r- 



E 

m,fc 



I'S'jo.fcl 



dwjo/djol iVj, 



(45) 



^ It is assumed hereafter that \du)j^ /djo | is not so small that the quadratic 
term in the Taylor expansion in I has to be retained. 



does not depend on the value of viscosity. This is physically rea- 
sonable because propagating disturbances decay before they reach 
at least one boundary. Accordingly, the calculation of the tidal re- 
sponse can be performed in a subdomain of the star with boundary 
conditions on the state variables corresponding to outgoing waves 
where these leave the subdomain. 

In Appendix B we support this conclusion by showing how to 
obtain J45b for a model of a non-rotating star, for which there is po- 
tential resonance with modes in a g mode spectrum, by a standard 
technique for solving the tidal response problem with appropriate 
radiation boundary conditions. 

It is important to note that the expression ( 1451 as well as the 
expressions ( 15 U below apply not only to the case of moderately 
large viscosity, but they also give mean energy and angular mo- 
mentum transfer rates obtained by averaging the factor Ak over the 
interval ^ 5 ^ 1. Indeed, from equation ( I44t it is easy to see 
that J dSAk = tv/k. Thus, these expressions are formally valid 
when orbital or stellar parameters change at uniform rates such that 
a large range of S is covered in the course of orbital and/or stellar 
evolution. 

Let us stress, however, that when the energy dissipation is very 
weak the orbital evolution may become stalled at some fixed value 
of 5 with the system located between resonances, see eg. Terquem 
etal(I998). 



3.5.2 The case of very low viscosity 

Now let us consider the case of small k <C 1. In this case the factor 
Ak ~ l/(«:^ + sin^ tt^/tt^), and we have from ( I42t 



E, 



^u.kJQl-^rn,k 



-47ry^ 

^ «fejo + [dujJdjoYsin^-KS/'K^ 



I'S'jo.fcl 



Ni, 



(46) 



When 5 is not very close to either zero or one the second term 
in the denominator of ( I46t is much larger than the first one. In 
this case the energy exchange is proportional to the value of vis- 
cosity as expected. However, close to resonances when uo^,kja > 
\dujjg/djo\\ sin7r5/7r| we can neglect the second term. In this case 
the energy transfer is strongly amplified being inversely propor- 
tional to LOv,kjo- 



3.6 Energy and angular momentum transfers expressed in 
terms of overlap integrals 

It is convenient to represent the quantities Sj^k entering into ex- 
pressions for the energy and angular momentum transfer and de- 
fined through equation ( I28t as an inner product of a quantity de- 
termined by the forcing potential and an eigenmode. Following the 
discussion above, the latter from now on will be taken to be those 
appropriate to the inviscid case v = Q for which we recall tj = ^ . 
Using equation ([8j for the Fourier components of the gradient 
of the forcing potential, we can write 



'j,fc 



27r 



d xpe 



-Qj, with Qj 



'™'^V(r^F2" 



imcp 



e, ■ V(rV2" 



(47) 



The quantities Qj , so defined, coincide with the so-called overlap 
integrals considered by e.g. Press & Teukolsky 1977, Ivanov & Pa- 
paloizou 2004, IP7 and others. Additionally, instead of using the 
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norm given by the expression ( 121) evaluated in the inviscid limit, it 
is more convenient to consider another norm 



n. = ^((^,l^,) + (^,|c^,)K 



(48) 



which reduces to a standard form in the limit of a non-rotating star. 
Let us express the expressions for the energy and angular momen- 
tum transfer in terms of Qj and Uj . Substituting ( 147 1 and l l48t into 
the expressions ( 135) . ( 139) and ( 140) we obtain 

7^ \ ^ \-^m,kQj\ 






Lc = 



\^m,kQj 



m — '-^ — — 

i^m.k Dk.j,j 
m,k,j 



and 



Ei ^ - 



l^y,kj{i^m,k + mQ) \Am,kQj 



m,k,j 



^m,k 



Dk, 



where 



(49) 



(50) 



does not depend on normalisation of eigen functions, and we use 
equation ( 125) to obtain the last equality. 



3.6.1 A dense spectrum of modes and moderately large viscosity 
In this case equation J45) together with ( 138) leads to 



E, 



-E- 



m,k 
E, = ^TT^ 






and 



m£l \ \Am,kQj 



tOm,k I \du}j„/dJQ\ 



(51) 



3.6.2 A dense spectrum of modes and very small viscosity 
In this case equation ( |46) similarly leads to 

Ec = ~ / , p''° l^m,fcQjo| , 



E 



<^'^,kJQ I 



IjJm.kD, 



and 



Lc — y '-ir^lAm^kQjol ! 



m,k 

E, = -J2i^ + 



(jJm.k I D, 



\Am,kQjo\ 



(52) 



m.k 

,2 I /J,. IJA \2„;„2„r/„2\ 



where D, = (^^,fcjo + (doj^a/djo) sin tt^/tt ). 



derived by IP7 Q As in the case discussed above they also can be 
represented as sums over quantities determined by the eigenmodes, 
which for a given m have eigenfrequency cok , that take the form 

a;|jSm(ajfc)P 



ABc = Stt^ ^ ■ 



Nk 



ALr. = 



AEi 



3 Y^ mUJk\Sm{i^k) 



m,,k 



N,, 



ujkjujk + mQ)\Sm{i^k) 

Nk 



(53) 



Representing again the quantities 5,71(0") as a product of two 
factors determined by the orbit and by the eigenmodes, respec- 
tively, we obtain a relation corresponding to (147) : Smio) — 
{l/2n)Am{(y)Qj. Using the normalised overlap integrals defined 
in l |50) we can rewrite ( 153) in a more convenient way as 



AEc = 2n''^\Arr^{uJk)Qk\\ 
m.k 

ALc = 2tv'^ y^ —\Am{iUk)Qk 



AEi 



2-K- 



n 



1 + !!!_') \Am{<jJk)Qk\^ 

^k / 



(54) 



Note that the quantities Am{o) entering ( 154) are discussed in e.g. 
Ivanov & Papaloizou (201 1) for the general case of a parabolic orbit 
inclined with respect to equatorial plane of a rotating staQ . 



4 ROTATIONALLY MODIFIED GRAVITY MODES 

4.1 Use of the traditional approximation 

In what follows we are going to apply the general theory discussed 
above to the situation when the energy and angular momentum ex- 
change between a rotating Sun-like star and its orbit occurs as a 
result of the excitation of internal rotationally modified gravity (g) 
modes. These are assumed to be associated with discreet low fre- 
quency spectrum. 

For the cases of a periodic orbit and of a nearly parabolic fiyby, 
the energy and angular momentum exchanges between the orbit 
and the star are characterised by two quantities associated with the 
noiTnal modes. These are their eigenfrequencies uij and the corre- 
sponding overlap integrals Qj. Thus for given orbital parameters 
and the rate of rotation of the star they determine the energy and 
angular momentum transfer rates in the case of 'moderately large 
viscosity' discussed above. In this regime the decay rate of a char- 
acteristically excited mode is small compared to its frequency as 
seen in the rotating frame. However, it is large compared to the 



3.7 Energy and angular momentum transfer in the case of 
nearly parabolic orbit 

The transfers of energy, AEc, AEi and angular momentum, ALc, 
between the orbital motion and noimal modes of the star during a 
periastron flyby on a formally parabolic orbit can also be expressed 
in terms of overlap integrals. For completeness we give expressions 
for these here. As mentioned above, in this case we should represent 
the results in terms of the Fourier transform of the tidal potential 
(see equation ^) instead of the Fourier series used in the previous 
sections. Expressions for AEc, AEi and ALc in terms of this were 



■^ Note that in equation (5) of IP7 there should be a (— ) sign in front of the 
expression on the right hand side and in equation (28) the imaginary unit 
should be inside the summation. 

■* Note that analogous expressions of Ivanov & Papaloizou 2011 differ by 
factor of two from those given in l l54t . This is due to the fact that in Ivanov 
& Papaloizou 2011, the summation over the azimuthal mode number is 
formally perfonned over positive and negative values of m, while in this 
paper we consider only positive azimuthal mode numbers. Note also that in 
this paper we take all quantities related to eigenmodes to be defined in the 
rotating frame, while in Ivanov & Papaloizou (201 1) the inertial frame was 
used. 
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frequency separation of neighbouring modes. Then the waves asso- 
ciated with the rotationally modified (g) modes are dissipated near 
a boundary, preventing the setting up of standing waves. This leads 
to a radiation boundary condition for the forced problem. However, 
these disturbances are solutions of the adiabatic linearised problem 
away from localised dissipative regions. When dissipative effects 
are very weak, an alternative prescription for their action on stellar 
perturbations should be used (see sections l3.5.2l and l3.6.2l 

The rotational angular frequency Q, can be expressed in units 
of the natural frequency of the primary star, fi* = yjGhh/R^, 
where Af, and i?, are the stellar mass and radius, respectively. As 
the non spherical distortion of the star is oc (f2/Sl,)^ it may be 
assumed to be spherically symmetric when (fi/fi, ) is small with 
internal structure the same as in the non-rotating state. 

The problem on hand has a symmetry property ujj — >■ —ujj, 
Q. — >■ — n. In the preceding sections it was assumed that Q, is posi- 
tive while the eigenfrequencies Uj can have either sign depending 
on direction of their propagation with respect to rotation of the star. 
In what follows it is convenient to adopt the alternative convention 
assuming that the eigenfrequencies are always positive while the 
case r2 > and Q < correspond to the modes propagating in the 
prograde and retrograde sense with respect to rotation of the star, 
respectively. 

When the orbital period is sufficiently large, tidal perturba- 
tions mainly interact modes of large radial order. In this case one 
can make use of the Cowling approximation which neglects vari- 
ations of self-gravity when describing free normal modes of the 
star and the 'traditional' approximation, where only components 
of the Coriolis force perpendicular to the radius vector from the 
centre of the star are included. When the traditional approximation 
is adopted and the star is assumed to be spherical the equations 
describing the stellar pulsation admit separation of variables (see 
below). The problem of finding ujj in such a setting has been con- 
sidered by many authors, either by numerical means or using the 
WKBJ technique ( e.g. Lee & Saio 1987, Bertholomieu et al 1978 
among others ). 



4.2 Calculation of the overlap integrals 

However, less attention has been payed to the calculation of the 
overlap integrals, Qj, Rocca (1987) attempted to calculate quan- 
tities called tidal resonant coefficients that are proportional to Qj, 
following the earlier work by Zahn (1970) who studied the same 
quantities in the case of non-rotating star 

However, significant potential inaccuracies remain, for exam- 
ple, Berthomieu et al (1978) assumed that the WKBJ approach is 
approximately valid in the convective regions of the star bounding 
a radiative region where the waves associated with the modes prop- 
agate. Although they did not consider the overlap integrals, this 
is not appropriate for the modes of low azimuthal order that are 
important for the tidal problem. In addition, when calculating the 
eigenfrequencies and eigenmodes, it has been generally assumed 
that, in a radiative zone, the square of the Brunt - Vaisala fre- 
quency behaves approximately linearly with the distance to the base 
of convection zone. Since for some models solar mass stars, this ap- 
proximation is not accurate (see Fig|4]below) we allow for a more 
general power law dependencqj 



Our approach to calculation of Qj is similar to the approach 
of Rocca (1987). However, in contrast to that work, we consider a 
general power law dependence, on the distance to the convection 
zone boundary, for the Brunt - Vaisala frequency. We revisit the 
calculation of ujj and Qj , mainly concentrating on calculation of 
Qj , since calculation of the eigenfrequencies is relatively straight- 
forward. We find two contributions to the overlap integrals, the first 
comes from the radiative region close to the base of convection 
zone and other comes from the convection zone itself. We correct 
some errors in earlier work, and find the asymptotic form at low 
forcing frequency. We compare the results we obtain with numeri- 
cal calculations for solar mass models in a regime where these can 
be performed. 



4.3 Linearised equations for normal modes 

Using the Cowling and traditional approximation in (|5), while ne- 
glecting dissipative effects and the external forcing potential, we 
get 



Lj^S, - 2it^ (n ■ f ) f X ^ + C^ = 



where 



Ci = -VP' 

p P^IP 



pi 

VP + iV^(^-f)f, 



p'(0 = -riPV-(^)-^-vp 

and the square of the Brunt - Vaisala frequency is given by 



AT^^i^f i^_ 



1 dP 



p dr \pdr TiP dr 



(55) 

(56) 
(57) 

(58) 



where Fi = {dlnP/d\np)ad, with the subscript ad indicating 
that the derivative is taken at constant entropy. 

We adopt spherical polar coordinates {r, 6, <j)) with ^r, S.e 
and ^0 being the components of the vector ^ in these coor- 
dinates, the component of fl acting in the radial direction is 
fl ■ r = (r2cosS,0, 0) where f is the unit vector in the ra- 
dial direction. It will be taken as read that all quantities of interest 
will be taken to be expressed in terms of Fourier series as specified 
as in (TJ with the equations applying to separate coefficients. Here 
p is the density and P is the pressure in the background state. The 
gravitational acceleration at radius r is g = GM{r)/r^ and M{r) 
is the mass enclosed within a sphere of radius r. 

The density and pressure perturbations are p and P respec- 
tively. We note that p , P and ^r satisfy the condition for adiabatic 
perturbations which gives 



2 ' 
CsP 



9 



ir 



(59) 



where Cs — y FiP/p is the adiabatic sound speed. Note that if 
A*'^ ^ everywhere as assumed in this paper, the operator C is 
positive and the linear eigenvalue problem can be brought into the 
self-adjoint form ( I17t . For simplicity we suppress the mode index 
(j ) below, wherever this doesn't lead to any ambiguity. 

As is well known (e.g. Unno et al 1989), both the forced re- 
sponse problem and the unforced (']/ = 0) problem considered in 
this section can be solved by adopting the separation of variables 
anzatz: 



^ Note that Provost & Berthomieu (1986) give WKBJ expre.ssion.s for the 
displacement vector in a radiative region valid for om' power law depen- 
dence. However, solutions in the convective are obtained numerically pre- 



cluding the development of an analytic asymptotic representation of the 
overlap integrals, which were not considered, at low forcing frequency. 
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^^ = S^— ' P 



together with 



p(r)H{e) 

2^ ' 



andP 



2-n 



(60) 



e{r) 



2;?(i- 



sine ^ ' dO I 



ie = 



i^ir-) fdH{9) 



2^ {1 ~ v^ cos"^ 9) \ de 



mvcoteH{e)\ 



(61) 



Here quantities with a hat are functions only of r, p = 20/ u 
and the factor Xj^phi is inserted for convenience (see below). We 
remark that the separation of the (^ dependence through the usual 
factor exp(im<^) is here taken as read. 

Substituting equations ( 160b and ( I61t into equations ( I55t -([59} 
we verify separability provided that the function Hid) is deter- 
mined through solving the eigenvalue problem specified by 

LH — —AH, where the operator L is defined by (62) 



d 
dfi 



(1 — v^jj,'^) dfi 



+ ■ 



mu{l + v fi 



(l-!/V^) V (l-i'V^) (1-M^) 



(63) 



with fj, = cos 9. The eigenvalues, A, specified by ( I62t . which de- 
fines the Laplace tidal equation, are also separation constants. The 
associated eigenfunctions are Hough functions ( see e.g. Longuet- 
Higgins, 1968, for a discussion of the eigenvalues and eigenfunc- 
tions). 

After separating the variables by the means as indicated above, 
it is found that the angular components of ( I55t require that 



C^ 



rcu^p 



(64) 



The radial component of ( I55t together with (I57t then lead to two 
ordinary differential equations for the radial dependence of ^(r) 
and P. These equations take the form 



d^ 
dr 

and 



'2 1 dP 



J^'^ p Uva cij 



dP , 2 



N^)i- 



P dP 



(65) 



(66) 



FiP dr ■ 

In general A is a function of v, and, accordingly, u. Therefore, 
when n 7^ the eigenvalues A and uj defined through the solution 
of equations ( |65t . ( 1661 ) and ( 1621 ) must be determined together. 

In the degenerate case Q. = Q, v = Q, ([62} reduces to Leg- 
endre's equation. Its solutions are accordingly are associated Leg- 
endre functions P™(/i) and A = /(/ + 1), where I is an natu- 
ral number. Since in this limit the tidal problem is such that only 
eigenmodes corresponding to Z = 2 are coupled to the dominant 
quadrupole component of the tidal potential, A can be set to 6. The 
eigenvalue problem for the determination of cj is accordingly de- 
coupled from the angular one in the non-rotating case. Equations 
(I65t-(I66|) then coincide with standard equations describing the free 
modes of pulsation of a non-rotating star under the Cowling ap- 
proximation, see e.g. Christensen-Dalsgaard (1998). 

4.4 Properties of the Laplace tidal equation 

We here summarize properties of equation (162) required for our 
puiposes. When the parameters p and m are fixed the operator L is 



self adjoint with respect to the natural inner product defined by inte- 
gration over (—1, 1) with respect to /i. There is a countable infinite 
set of solutions coiTesponding to different characteristic values of 
A, Ai. The corresponding eigenfunctions. Hi, are orthogonal and 
normalised such that 



dfiHi{fi)Hk{p) = Sik. 



(67) 



In general, the eigenvalues , Ai , can be either positive or neg- 
ative. In this paper we assume below that they are non-negative, 
which is necessary for the associated solutions to correspond to 
wave-like perturbations in stably stratified regions, as required for 
potentially resonant rotationally modified gravity modes. Thus the 
Ai, i — 1, 2..., will be arranged such that A^+i > A,; with Ai cor- 
responding to the minimum non-negative eigenvalue for a given u 
and 771. 

As stated above, when i/ = 0, the Hi are proportional to as- 
sociated Legendre functions, P;'". In the general case, f ^ 0, the 
functions Hi may be represented as an infinite series of associated 
Legendre functions through 



H,. 



n^O 



"i,nP('?n + 2), with P(™„ + 2) = ^J^i."-'^^- (68) 

i>0 



Here the scaled associated Legendre functions 



prif^) 



{2l + l){l~my. 



2{l + m) 



pr{p) 



(69) 



satisfy the condition J dfiP[^{fi)Pl!^{fj,) — Si^ii. The second 
equality in {68} follows from orthogonality and normalisation of 



each of the sets of functions Hi and PJl 



. We calculate the de- 



(2n + 2)- 

composition of the tidal potential as a series of Hough functions 
using equation ([68} truncating the summation at sufficiently large 
value of, n, with help of the procedure described in Ogilvie & Lin 
(2004). 

As the dominant tidal potential has the angular dependence of 
a quadrupole, it follows from equation ( I47t . that the overlap inte- 
grals involve a product of P™ and Hi, and therefore, the coeffi- 
cients 



ai{v) = OLifi 



dfiH.P? 



(70) 



play an important role. It describes the 'angular' coupling of the 
rotationally modified g-modes with the tidal field within the frame- 
work of the traditional approximation. 

For moderate values of i^ ( for which, lu] ^ 3, say) only 
the quantities cortesponding to i = 1 are important. We calcu- 
late Ai (u) and ai (z/) numerically as well as an additional quantity, 
Ie{u), determined by the Hough functions and defined below (see 
equation ([84} and plot them in Fig.[T]for m = 2 and in Fig.[2]for 
m=0. Note that in the latter case from the definition of Ai {v) and 
Qi iy) it follows that they are even functions of v while Ig is an odd 
one. Therefore, we show only positive values of v in Fig. [2] 

The case m = 2 is the most important for tidal interaction. 
We see from Fig.[T]that Ai (v) and —Ig (v) increase as i/ decreases. 
On the other hand a\{u) increases with v until v = —1, levelling 
off at larger values. As a\ (v) measures the strength of coupling to 
the tidal field, this latter behaviour results in the stronger excitation 
of modes associated with positive values of v. 

The values of v are in general different for different modes 
since these are obtained from the simultaneous solution of eigen- 
value problems in the radial and 9 directions. Nonetheless, we clas- 
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Figure 1. The values of Ai (solid line), ai (dashed line) and —Ig, where Ig 
is defined through equation l l84t (dot-dashed line) for m = 2 as functions 
of the parameter u. 



sify the eigenmodes using two indices, the index (i) coiTesponding 
to a particular A; {i^i), which accordingly indicates the order of the 
angular eigenvalue, and the index (n) describing the number of 
nodes in the radial direction, for a given m. Thus our mode index 
j, used in previous sections, should be understood as a pair (i, n). 
Note that when the star is non-rotating, and, accordingly, 
u = 0, from equationdM} it follows that ^ = ^Y^^Br + ^^rWr, 
where e,- is the unit vector in the radial direction. 



4.5 An expression for the overlap Qj in the low frequency 
limit in the traditional approximation 

We now reduce the expression ( 147) for the overlap integral Qj to 
a form involving an integral over the radial coordinate only in the 
traditional approximation. This is useful for either performing nu- 
merical calculations or undertaking further asymptotic analysis (see 
below). We work on the expression for Qj given by equation l l47t . 
It is convenient to express the term V{r'^Yp) that appears 
there as 2rY^er + r^VY^T ■ We then split the integral for Q-, into 
two parts, the first containing 2rY'J^^* ■ r and the second contain- 
ing r'^^'j ■ \7Y^. We remark that the forms of both the free and 
forced problems specified by (155) are such that after separation of 
variables, the functions ^ and ^ written as ^j and ^j, when these 
correspond to a free normal mode j, may be taken to be real. Not- 
ing that Y^ = P^ exp(im<;/))/^27r and using {TOll, the first part 
of {47} can be transformed to become 



2ai 



r drp^j. 



(71) 



To deal with the second part we first note that equation ( 163) is 
equivalent to the relation 

rVx(e''"'*^,) = -A^^fH^e'""^/V2^, (72) 

where Vx = V — ff ■ V. Making use of ( 172) . performing an 



Figure 2. Same as Fig. \T\but the case ; 
line shows Ig . 



is shown and the dot-dashed 



integration by parts, and noting ( 168) and ( 170) . the second part can 
be reduced to the form 



ctiAi 



Rt 



3 J ^S 

r drpt^j . 



Summing the contributions of the two parts, we obtain 



Qj = aiQrj,v7ith Qr 



r^drp{2i,+Kil). 



(1'i) 



(74) 



When the star is non-rotating q; = 1, A = 6, and the expression 
( 174) reduces to its standard form, see e.g. PT, IP4. 

We remark that, when considering which regions of the star 
contribute to the overlap integrals, the original expression l |47) can 
also, after a straightforward integration by parts, be reduced to an 
integral for which the integrand is proportional to the mode density 
perturbation. Viewed in this way, only regions with non vanishing 
density perturbation contribute, consistent with expectation for a 
tidal interaction controlled by the gravitational potential due to a 
perturber. 

The expression for Qrj can be transformed into a form that is 
explicitly oc uj^ which is thus very convenient in the low frequency 
limit, see also Rocca (1987). In order to proceed we first use (164) 
to express P in terms of ^ and substitute the result in equations 
( 165) and | |66) . These then yield the two relations: 



1^7 / rN q , c. 






(75) 
(76) 



where a prime indicates differentiation with respect to r. Now we 
substitute the right hand side of eauation( l76) for ^j in ( 174) and then 
integrate the teiTn proportional to [r^ p£,)j by parts, assuming that 
^j is smooth, but we allow for the density p to undergo a discontin- 
uous jump at some point r = r, as is found to occur in some stellar 
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models. In addition we assume that there are no boundary contribu- 
tions coming from the surface or centre of the star. We accordingly 
obtain 



Q 



Rt 



r drp 






(77) 



(«.|C^.) = ^ 



\ 9 

where the term in square brackets in ( 177) should be understood 
as being equal to the difference between the expression in square 
brackets taken at radii r, + e and r* — e and then taking the limit 

In order to get an expression for Qr that is explicitly propor- 
tional to a;| we substitute the right hand side of equation fI5l for 
^j into l |77| l and then integrate the term proportional to {r£,j) by (^j|B^j) = —Qlrle 

parts assuming the same conditions as above. We then get 



Using equation J56b to specify J81b in terms of the linear pertur- 
bations, making use of I l57t -(l66b. perfoiTning an integration over 
the angular coordinates by parts, and then making use of ill) , we 
obtain 



r^drpi^^+M^ff). 



(82) 



This expression reduces to the standard one in the non-rotating case 
when A = 6. Using l ll2t . with ft replaced by its radial component 
{ft ■ f )f as appropriate for the traditional approximation, together 
with l l61b . we obtain 



(83) 




where the surface term 



S^ = r% [p]^ 



2r 



-3 g [Pcf ] 



Sl± 



(78) 



(79) 



In order to evaluate S^ we take equation l l75t , multiply both sides 
by pTV^ and substitute the expression given by ( I58t for the square 
of the Brunt - Vaisala frequency. Integrating the result between 
r, + e and r, — e, making use of the fact that the displacements are 
are finite at r = r* and that the pressure gradient is equal to —pg, 
we readily obtain S — 0. 

So far the expressions for Qrj, though adapted to low fre- 
quencies, are exact and can be used for numerical evaluation. To 
facilitate analytic discussion in the relevant limit of low frequency, 
we express the integrand in ( 177) in terms of ^j only, and in so doing, 
obtain a simplified expression for Qrj that neglect teiTns of order 
u^ and higher. 

In order to do this, we substitute the right hand side of equation 
([76} into l l78b . first noting that the last term may be neglected to 
the order we are working. Integrating again by parts, assuming no 
boundary contributions or internal jumps in ^® , we obtain 

cR, 
LOi I r drph{r)^j + 0{u]j), where 



h{T) 



1/ (12 ^4i;) 
5 V ^« 



1 r^g 2{rgY , {&-v)rg 



A. 



S2 



(80) 



4.6 An expression for the norm 

In order to calculate energy and angular momentum exchange rates 
resulting from tidal interaction, in addition to the overlap integrals 
Qj, we also require the norms n-, ( see eg. section ITSt . Here we 
develop a form for Uj in terms of one dimensional angular and 
radial integrals that can be readily evaluated. To do this we trans- 



form the expression ( |48b for the norm, Uj , so that it involves in- 
tegrals only over the radial coordinate. Using the that ^ satisfies, 
C^ + uu B^j , ( 148b can also be written as 



2/- 

w 4 



27r 



= ^((^.|c^. 



+ f(^.ie^.))- 



(81) 



where 



Ir = I r^drp{^j)'^, and le 



fidfi 



up{l - p ) 
Substituting 



dH 



dp 



+ m{l + u p )H-— + — — 

dp L — p^ 



(84) 



and l l83b in J81b we get rij — Ust + rir, 

Rt 



where Ust = I r drp{£,j + Ai(^j ) ), and Ur = vlrh- (85) 





We remark that the integrals with respect to r can be straightfor- 
wardly evaluated once the eigenfunctions are known. The angular 
integral Ig is a function of v only. The dependence of —Ig onv is 
illustrated in Fig.[T] 



ROTATIONALLY MODIFIED GRAVITY MODES FOR 
MODELS OF SOLAR-TYPE STARS 

In this section we describe the numerical evaluation of the eigen- 
frequencies u}j and the overlap integrals Qj for rotationally mod- 
ified g -modes under the traditional approximation for stars with 
solar-type structure. We also highlight aspects of the models that 
are relevant to a subsequent asymptotic analysis. 



5.1 Stellar models considered 

Our calculations are for a rotating star with a radiative core and 
extended convective envelope. The structure is assumed not to be 
modified by rotation, such modifications being a second order ef- 
fect. This has the advantage that standard spherical models can be 
used. 

We have considered an almost zero age model of the star with 
M* = Mq Ri 2 • lO^^g and age t ^ 1.66 • 10** yr provided to 
us by I.W. Roxburgh (see Roxburgh 2008 for a description of his 
numerical code) referred hereafter to as IR model, and also a model 
of the present day Sun (Christensen-Dalsgaard et al. 1996) with age 
t = 4.6 ■ lO^j/r, referred hereafter to as the CD model. Plots the 
density p and the square of the Brunt - Vaisala frequency, A^^, 
against radius, for these models are given in Fig. |3] Note that the 
total stellar radius of the IR model, R, is smaller than that of the 
present day Sun, taken to be _R, = Rq = 7 ■ 10^° cm. For the IR 
model we have 7?, = 6.3 • lO^^cm approximately. 
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Figure 3. The density p (in units of the mean density 3M*/(47rR2) and 
square of the Brant - Vaisalii frequency, A'^^ (in units of GMt/Rt) as 
functions of the radius r expressed in units of R« . The solid curves corre- 
spond to the CD model while the dashed ones are for the IR model. The 
curves monotonically decreasing with r are for the density distributions, 
while the curves having maxima at some values of r are for the Brant - 
Viiisala frequencies. 



5.2 Numerical approach 

The eigenfrequencies are obtained by finding appropriate solutions 
of ( I65t and l l66t . in the core and in the envelope and applying 
matching conditions at the base of convection zone. Since we con- 
sider low frequency modes we can neglect the last term in the 
brackets multiplying P in ( I66t as this is important only very close 
to the surface in the limit o;^ — >■ of interest. The radial eigenvalue 
problem is solved together with the azimuthal eigenvalue prob- 
lem. To solve the former the standard fourth order Runge-Kutta 
method was used. To solve the latter for a given eigenfrequency, 
we adopted the matrix method of Ogilvie and Lin 2004. We em- 
ployed ~ 10® — 10^ grid points in the radial direction in order to 
ensure numerical convergence for eigenfrequencies uj <~ 0.1 in 
natural units. Since the CD and IR models were provided with a 
much smaller number of grid points, 1282 and 2001 respectively, 
linear interpolation was employed to supply the data for our finer 
grid. Having obtained the eigenfunctions, we then use the expres- 
sions ( I78t and ( I84t to calculate the overlap integrals and norms. 



5.3 Asymptotic analysis in the limit cj 



0. 



In subsequent sections we shall discuss the asymptotic form of the 
normal modes, overlap integrals and norms in the low frequency 
limit following a WKBJ approach. In order to proceed we first dis- 
cuss the behaviour of A'^'^ in the vicinity of the interface between 
the radiative and convective regions of the stellar model, this being 
important consideration for the asymptotic analysis . 



Model 


A 


q 


Tc 


CDl 


4L54 


0.747 


0.722 


CD2 


24.32 


0.617 


0.721 


IRl 


41.4 


0.737 


0.731 


IR2 


23.1 


0.564 


0.73 


CDL 


106 


1 


0.724 


IRL 


110 


1 


0.733 



Table 1. Results of fitting of the dependence of N"^ on r for the CD and IR 
models. Different entries in the table correspond to either general fits made 
with different data sets, or fits assuming q = 1 , as described in the text 
(see also Fig.|4). Note that A and r^ are given in natural units such that r 
and A^^ are expressed in the units of i?» and the natural stellar frequency 
GMt/Ri, respectively. 



5.4 Behaviour of the Brunt - Vaisala frequency near the base 
of the convection zone 

In previous studies (e.g. Zahn 1970, Berthomieu et al 1978, Rocca 
1987) it was assumed that near the boundary between convective 
and radiative zones, but on the radiative side, 7V^ oc the difference 
X = \rc — r\ ^ r, where Tc is the radius of the boundary , defined 
by the condition N'^{rc) = 0. However, this is not to be the case 
for many stellar models (see Barker 2011). In particular, as seen 
from Fig.|4] for the solar models we use, N'^ decreases more slowly 
with X than expected for a linear dependence. Here we consider 
more general possibilities, by assuming that when Vc > r, so that 
X = Tc — r > 0, for a; ^ r we may write 



N^ 



Ax" 



(86) 



where g is a constant, see also Provost & Berthomieu (1984). Note 
that as is seen from Fig.|4]the dependence of A'^^ on x = r-c — r is 
not strictly given by a single power law. Local fits result in smaller 
values of q at larger values of x. However, this feature is not very 
important for our purposes. For our estimates of the overlap inte- 
grals below, only the local value of q that is found to be applicable 
to the range of x such that 1 > A*'^ > 0.05 is important. This 
is needed to find the contribution to Q,. given by equation l lll2t . 
which plays a role only when eigenfrequencies are in the range 
cj ~0.2- 1. see Figs. |7l9l 

We use equation l |86| l to fit the dependence of the Brunt - 
Vaisala frequency on r in the CD and IR models in the vicinity 
of the base of convection zone using the Levenberg-Marquardt al- 
gorithm. 

Since the dependence is not, in fact, an exact power law we made 
fits using two different data sets in order to explore the range of 
the resulting values of A and q. For the first data set the number of 
data points is constrained approximately by the condition A*''^ < 2, 
while for the second, and more extended data set, the number of 
points is constrained by the condition 7V^ < 4. Additionally, we 
made fits to the data fixing the power index g = 1 for the data sets 
with A'^^ < 2. The results are illustrated in Fig.|4]and tabulated in 
tabled 

The results labelled CDl, CD2 and CDL are respectively for 
the fits using the data sets with A'^^ < 2, A'^^ < 4, and with assumed 
linear dependence of A'^^ on x, for the CD model of the star The 
labels IRl, IR2 and IRL describe the same type of fits made for the 
IR model. 

From Fig.|4]we see that the fits obtained with a smaller number 
of data points give somewhat larger value of q « 0.74 whereas the 
larger data sets are fitted quite well by ( 186) with q ~ 0.6. This 
indicates that the value q required for the most accurate fit may 
decrease with the radial extent of the region of interest. This in turn 
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Figure 4. N'^ close to the base of convection zone plotted as a function of 
radius, r, expressed in units of ij* . The filled circles show the data points 
used to make the fits. Those with larger values of A^ ^ at a given r correspond 
to the IR model. The solid and dashed curves show the CDl, IRl, and the 
CD2, IR2 fits respectively. The dashed straight lines show the fits CDL and 
IRL which had q = 1, see table[T]for the values of the parameters A and q 
appropriate to these fits. 



may be a function of frequency when normal modes are considered. 
We discuss below how quantities obtained in our analysis depend 
on the assumed fit. 



6 CALCULATION OF EIGENFREQUENCIES, 

EIGENFUNCTIONS, OVERLAP INTEGRALS AND 
NORMS FOR HIGH ORDER ROTATIONALLY 
MODIFIED GRAVITY MODES FOR MODELS OF 
SOLAR-TYPE STARS FOLLOWING A WKBJ 
APPROACH 

In order to provide a check on numerical calculations and to extend 
results to frequencies below which the numerical approach can run 
into difficulties on account of inadequate resolution, we perform an 
an asymptotic analysis following a WKBJ approach. We first obtain 
expressions for eigenfrequencies and eigenmodes and then go on to 
use them to evaluate overlap integrals and norms. 

6.1 Asymptotic solution of the radial eigenvalue problem 
providing expressions for the radial eigenfunctions and 
associated eigenfrequencies 

We now discuss the asymptotic form of the eigenfrequencies and 
eigenfunctions found from solving the radial eigenvalue problem 
in the low frequency limit adopting a WKBJ approach. In what 
follows, for simplicity we make the assumption that Fi — const 
throughout the regions of the star of interest, setting Fi = 5/3 in 
numerical expressions. Note that these assumptions can be easily 
relaxed if needed. Similarly the term involving l/cj in ( 165b be- 
comes important only very close to the surface it be neglected in 



the same limit. This term can also be taken into account in a more 
advanced variant of our scheme. 

6.1.1 The WKBJ solution in the radiative core 

Under the assumptions mentioned above equations d65t and J66b 
can be combined to yield the second order differential equation 



An approximate WKBJ solution in the radiative core, sufficiently 
far from the centre that N'^ /ijp — 1 may be replaced by N'^ /uj'^ in 
l l87b . has the standard form (see e.g. Christensen-Dalsgaard 1998) 

C = CwKBjr~''^^^{pNy^''^ s\-a{W/u) + (t)wKB.j), y^^are 



W 



A 



1/2 



dr 



N, 



(88) 



CwKBj is a constant and the angle (Jjwkbj is fixed by the condi- 
tion that the solution matches onto a regular one as r — >■ as 



In 



bwKBj = ---, with I = -(VT+IA- 1), 



(89) 



see e.g. Vandakurov (1968jj. Note that I is not, in general, an in- 
teger When the star is non-rotating we recall that A = 6, and, 
accordingly, 1 = 2. 

6.1.2 The solution close to the base of convective envelope 

We look for a solution that is valid in a radial region of extent r <^ 
Tc in the radiative region close to its boundary with the convection 
envelope. To do this we substitute J86t in J87t . neglecting unity in 
the brackets in the last term, the variation of P and p, and recalling 
that X = Vc ~ r. Equation dSVt then reduces to 

d^^ AAx\ „ 
dx^ riiij-' 

The general solution of ( I90t is a linear combination of Bessel func- 
tion Jv[z) such that 



i 



s/x(C^J 1 iz) + C7J 1 iz')). where 

g + 2 g + 2 



K 3+1 
— X 2 , 



K = 



(91) 



(g + 2)rc' 

Note that the expression ( I91t is equivalent to what is obtained in 
Provost & Berthomieu (1984) in the limit of small x provided that 
one discards higher order terms as discussed in their paper 

The solutions l |88b and l |9U must be matched in a region where 
x/r <C 1 but z ^ 1. This is easily accomplished by the stan- 
dard method of comparing the former with the latter after utilising 
the asymptotic expansion of the Bessel functions for large z. The 
matching procedure allows us to express the constants Ci and C2 
in l |91t in terms of Cwkb.j through 



Ci = 
C2 = 



nK cos 



■kK cos( 



{ 4>c-nq/{A{2 + q))) 
sin(^/(2 + g)) 

.-^(4 + g)/(4(2 + g))) 
sin(^/(2 + g)) 



(92) 



^ It can be shown that close to the centre we have approximately: £,core. oc 
r""^' ^ J;^i ;2(Ccore^), where Ju{z) is a Bessel function and Ccore is a 
constant. This solution is regular when I ^ 1. At sufficiently large r the 
solution should be matched to |88j. This determines the phase 4>WKBJ- 
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where C, = CwkbjPc r^ A ^''*, and 



= (pWKBJ + 



%/A 



■= dr 
— N. 
r 



(93) 



But note that the approximations that lead to equation l |90t 
imply that this equation must be modified when A^*^ ^ uj^, which 
corresponds to 



X < Xc 



A 1UJ1 . 



(94) 



On the other hand from ( |91t it follows that when z <^ 1 the solution 
approximately satisfies d?^/dx^ = and the neglected term can- 
not produce significant changes over the small shrinking domain 

{Q,Xarit) aSW — >• 0. 

6.1.3 The solution in the convective envelope 

Extending the solution into the convective region, we can set ap- 
proximately A'^^ — 0. The stellar gas, therefore, can be approxi- 
mated as having barotropic equation of state with P — Kcp^^, 
which corresponds to a polytrope with polytropic index equal to 
1.5. Assuming that the convective envelope has a negligible mass 
compared to that of the star, of mass A/* , the hydrostatic balance 
equation takes the form 

dP _ GM,p 
dr r'^ 



(95) 



Solving l |95t to find the dependencies of the pressure and density 
on r, assuming these vanish at r = _R, , we obtain 



P = P* 



R, 



- 1 



with p, = I %:i)g^ I .(96) 



Setting N'^ = in equation ([87}, using the polytropic equation of 
state and equation l |96t we obtain 



1 



d$e 



dr^ 



(Fi- l)r(l-r/i?*) dr 



u. 



0, 



(97) 



where <I>c = r^pS,- We look for a solution to ( |97t , which tends to 
zero when r ^ Ri,. This solution is expressed through the hyper- 
geometric function, F{a, /3, 7, z), as 



^c = Ccr^x^i-^ F{bi,b2,b3,x), 
where x — 1 — r/R^, 



bi,2 
63 



(98) 




2ri -1 



ri 



and 



and C'c is an arbitrary constant. As we discuss in the next section, 
for our purposes it is important to know the quantity 



(100) 



calculated for the solution discussed above, at the base of the con- 
vective envelope, as a function of A. This quantity can be calculated 
from ( |98t . It will be seen from the above analysis that this quantity 
depends on the eigenfrequency only through its dependence on A. 
The dependence Be on A, for Fi = 5/3 and rd R* = 0.7, is 
illustrated in Fig.|5] 




Figure 5. The quantities Be JlQOt and 10^ i"c (see (TTS)) determined from 
the solution j98t as functions of A for r^ = 0.7 R*. Be is shown as the 
solid curve. Note that the values of Be and le found from the asymptotic 
analysis are in quite good agreement with those obtained using numerically 
calculated eigenfunctions. 



6.1.4 Matching conditions and eigenfrequencies 

In order to mach the solutions in the radiative zone and the con- 
vective envelope we use the solution l |91| l in the limit 2; — )■ to 
calculate —{d^/dr)/^ as the base of the convective envelope is 
approached from below. Then, we equate the resulting expression 
to ilOOt which gives the same quantity as the base of the convec- 
tive envelope is approached from above. This gives a compatibility 
relation, which will be satisfied only when u) is an allowed eigen- 
frequency. In the limit 2 — >■ equation l |91t gives 



i 



Ci{K/{2uj)y^^'+''^ 



C2 {K/{2cu)) 



-l/(2 + q) 



r{{3 + q)/{2 + q)) 
1 



+ 



F((l + g)/(2 + g))' 



(101) 



where F(2:) denotes the gamma function. We then find 
~{d£^/dr)/^ in the limit z —¥ Q from dlOU and equate it to the 
same quantity determined from the solution in the convective enve- 
lope given by equation dlOOb . 

Equation l |92t can be used to specify the ratio C1/C2 which, 
after its elimination from the expression, yields a compatibility con- 
dition from which the eigenfrequencies can be determined. This 
compatibility condition can be written down as 



3(<^e - 7rQ/(8 + 4g)) 



F((3 + g)/(2 + g))^,_ . 



-^(4 + g)/(8 + 4g) F((l + g)/(2 + g)) 



Bcw5+i,(102) 



where uj — 2u)/K, and 4>c is given by equation l l93t . 

In the limit oj — >■ the term on the right hand side of ( 1102) 
is small. This allows us to solve equation \W1\ iteratively. First we 
set the term on the right hand side to zero. [j. Then we find that we 



This is equivalent to setting Ci = which con'esponds to the radial 
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Figure 6. The difference 5 = \(^WKBJ ~ ^n\l^n between the WKBJ 
eigenfrequencies it^wKB.] and the corresponding frequencies determined 
numerically tj^. The results are plotted as a function of ojn calculated 
for the CD model of the present day Sun and expressed in units of 
yjGMf/R^. The value of / defined in equation 1104) is approximately 

equal to 14.91^ / „3* ■ Since all the fitting models shown in table[T]give 

quite similar curves we show only the results corresponding to the CD2 
model. See the text for the description of different curves. 

must have (f)c = 7r(n+ l/2 + q/(4(2 + g))), where n is an integer. 
Using this result in equations ( 189) and ( 193) we obtain 

•/Xl 



'■ OJo 



7r(n + (1 + 0/2 + 5/ (4(2 + g)))' 



where 



dr 



N. 



(103) 



(104) 







Now let us take account of a small non zero Be by setting 
<j)c = 7r(n + 1/2 + 9/4(2 + q)) + A<^, where Aijf) is assumed 
to be small. We substitute this expression into ( I102t and solve the 
resulting equation for A(^ using perturbation theory to linear order, 
thus obtaining 



and accordingly 



r((3 + g)/(2 + g)) 
r((l + 5)/(2 + g)) 



^^KI 



7r(n + (1 + 0/2 + ?/(4(2 + g))) + A0(a;o) 



i(7r/(2 + g)),(105) 



(106) 



We recall that for rotating stars equation (1106 1 gives only an im- 
plicit equation for eigenfrequencies since in that case A is a func- 
tion of CJ. 

We remark that for low frequencies the influence of the con- 
vective envelope solution on the spectrum, through Be, provides 



eigenfunction having its derivative close to zero at the base of convection 
zone 



only a small correction to the eigenspectrum. This is consistent 
with the idea that the properties of the rotationally modified g mode 
spectrum are almost completely determined by the radiative zone. 

We now compare results obtained from our analytic expres- 
sion for the eigenfrequencies given by equation ( 1106) with those 
obtained from an expression for ui given in Berthomieu et al (1978) 
and also with those obtained from numerical calculations for the 
CD model of the Sun in Fig.|6] The simple expression obtained by 
Berthomieu et al (1978) differs from equation ( 1106) by replacing 
the denominator in that equation with 7r(7i+ i(l + ^/A + j)). 

In Fig. [6] we show the quantity 5 — \ijJwkb.j — cj,i|/cj„, 
where ijJwkbj and ujn stand for the eigenfrequencies determined 
from the WKBJ analysis and numerical calculation respectively, as 
a function of lj„ expressed in units of yjGMt/R^. The curves 
with symbols showing positions of the eigenfrequencies are calcu- 
lated from our analysis, while the curves with no symbols are for 
eigenfrequencies determined from the expression of Berthomieu et 
al (1978). The solid, dashed and dot dashed curves are for SI = 0, 
n = 0.420* and Q. — — 0.42S1, , respectively. It is seen that the de- 
viation between the numerical and analytical results does not vary 
monotonically with S increasing as to becomes very small. This 
feature is most probably a consequence of numerical errors. The 
curves corresponding to fl = and Q — 0.4251, are similar to 
each other, this can be explained by the rather weak dependence of 
A on ly when ly = — is positive, see Fig. [T] Since the WKBJ ap- 
proximation gets better with increasing numbers of nodes n, which 
is, in its turns, gets larger with increasing A, the deviation 5 cor- 
responding to those curves is appreciably larger than the devia- 
tion corresponding to retrograde rotation at w < 0.8 because A 
increases significantly with decreasing uj in that case, see Fig.[T] 

As seen from Fig. |6] our approach gives deviations from the 
numerical results that are smaller than those of Berthomieu et al 
(1978) for prograde rotation and no rotation. However, the differ- 
ence is typically rather moderate being of the order of several tens 
per cent. For the case of retrograde rotation and large values of 
ij the Berthomieu et al. (1978) expression gives somewhat better 
agreement though the difference is small. When uj < 0.847f7, our 
approach is closer to the numerical one than that of Berthomieu 
et al. (1978). Note a non-monotonic behaviour of the curve cor- 
responding to the Berthomieu et al (1978) expression close to 
UJ = 0.847f2,. This value of uj corresponds to i^ = 1, which is 
a special value of this parameter delimiting the modes dominated 
by the gravity (f < 1) and Coriolis forces (f > 1). We suspect 
that this feature may be purely numerical because of the behaviour 
of the singularities of equation ( 162) . determining A, as v passes 
through unity. 

For our estimates of the overlap integrals given below it is 
convenient to substitute the explicit expression for the phase (j)c 
calculated above into equations ( 192) for the constants Ci and C2. 
Making use of ( 1105) to do this, we easily obtain 



,2^x2/(2+9) r((3 + g)/(2 + g)) 



^' - (-') V S^^ (^) r(l + g)/(2 + g)) 



o* , 




C2 = (-l 

r((3 + g)/(2 + g)) 
r((l + g)/(2 + g)) 



-^(i^Bj'-^r'''\ 



s(n/{2 + g)) a 



(107) 
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6.2 Asymptotic expressions for the normalised overlap 
integrals Qj 

Having determined eigenfrequencies and eigenfunctions, we can 
use them to estimate tlie noiTnalised overlap integrals determining 
the strength of tidal interactions. We evaluate the unnormalised Q 
derived from ( 147) . and the norm n given by l |85K from now on we 
omit the mode subscripts {i,j) below unless absolutely necessary). 

The overlap integral is written as Q = aQ,.. The factor a, 
coming from the angular integration is calculated numerically and, 
as it does not depend on the radial foiTn of the eigenfunctions, it is 
assumed to be known function of frequency (see Fig.[T). As we are 
working in the low frequency limit, we use the reduced form ( |80t to 
evaluate the factor Qr- It is important to note that the evaluation of 
this expression can be simplified for solar type stars. For these, the 
main contribution to the integral comes from, either the transition 
region between the radiative core and the convective envelope, or 
from the convective envelope itself. We discuss these contributions 
separately below. Note that in view of various integrations by parts, 
these contributions do not necessarily relate to physical contribu- 
tions from the regions concerned. 

In what follows we adopt natural units, expressing radii in 
units of the stellar radius R,, mass in units of Af, , frequencies 
in units of y/GMt/R, and density in units of the stellar mean 
density p* = 3/(47r){M«/J?J). In these natural units the factor 
CM* in dlllb can be set equal to 1. The normalised overlap inte- 
grals have dimension of y/MZR, . Accordingly we express them in 
these units. From equation ( I47t . given the fact that we express the 
density in units of the mean stellar density, it follows that expres- 
sions for the normalised overlap integrals obtained with quantities 
expressed in our units, must be additionally multiplied by factor of 
•y/3/(47r). This will be assumed implicitly. 

6.3 An expression for the norm 

We consider the expression ( 185) which gives the norm n written 
as the sum of two contributions such that n = rist + n^. In the 
limit Lj — )■ 0, the main contribution to the radial integrals in (185) 
is expected to come from the region in the radiative core where 
the radial wavelength and group velocity are smallest and accord- 
ingly radial displacement is given by the WKBJ expression l l88b . 
In the same limit and in the same region of the star, |^ | 3> \^\ as 
expected for low frequency modes in a stratified medium. Thus it 
follows from equations ( 164) . ( 165) and ( 188) that 

CwKB.J 



e 



'^WKB.J 



Al/2^ 

W 



-p'^^^r-^/^N^^^ cosi'iwKBj) with 



(108) 



Substituting this expression into the integrand determining Ust in 
( 185) , neglecting the contribution proportional to ^ there and adopt- 



ing the average value of cos {"^WKBj) 



we get 



Ust 



-.2 
-■WKBJ 



2lj2 



(109) 



where / is given by equation ( 1104) . Since n^ — ulrh, with the 
integral Jr, defined through ( 184) , readily being seen to be approxi- 
mately equal to rist/A on account of l^'^l ^ |^|, we finally obtain 



Ust + rir 



^WKBJ r / 1 I '^-'9 

2a;2 ' ^ 



T 



(110) 



where le is defined through equation ( 184) . We recall that in general 
the latter has to be calculated numerically. 



6.4 The estimate of the overlap integral and Qr 

In order to calculate the integral in the expression for Qr, we use 
the approximate expression ( 180) . The integrand in that expression 
is proportional to the function h{r) defined there. We remark that 
the main contribution to the integral comes from regions close to 
the boundary between the radiative zone and the convective en- 
velope where the eigenmodes have their largest scale. The term 
proportional to the Brunt - Vaisala frequency in ( 180) is small in 
these regions and we set in to zero. Also, as indicated above, the 
mass of the star may be taken to be constant in these regions so 
that g ~ GMt/r^. Substituting this into the expression for h{r) in 
( 180) , we obtain the useful and simplified form 



h{r) 



GM. 



:('-?) 



(Ill) 



The integral may be regarded as being composed of two con- 
tributions, which should be evaluated separately, the first comes 
from the radiative region close to the base of convective envelope, 
and the second from the envelope itself. 



6.4.1 The contribution to the integral from the radiative region 
close to the base of the convective envelope 

In order to calculate this contribution we set all quantities in the in- 
tegrand, multiplying the radial displacement ^, equal to their values 
air = Tc and assign the subscript (c) to them, e.g. pc = p{rc). We 
assume that the distance from the base, x — re ~ r > 0, is small 
and change the integration variable from x to z as defined in equa- 
tion J91b . We substitute ^ as given by ( 191) into the integral in ( 180) . 
Taking the lower limit to be zero and formally extending the upper 
limit of integration to infinity gives the contribution from the radia- 
tive region. Proceeding in this way we get a linear combination of 
two integrals containing Bessel functions and powers of z. These 
integrals can be evaluated by standard methods using Gradshteyn 
& Ryzhik (2007) giving the contribution to Qr from the radiative 



zone, Qrr as 



= 1- 



30 \ 23/(^+^) .3/2 



c 



r(2/(2 + g)) 



C2 



r(i/(2 + g)) 



'r{(l + q)/i2 + q))^ V(g/(2 + g)) 



(112) 



where Dc = (cj/iC)^/<^+«^ and we recall that K is defined in ( |93l l. 
By construction this expression does not include the contribution 
to the integral from the convective envelope which is considered 
separately below. Next we substitute Ci and C2 as given by ( 1107) 
and use ( 192) to replace C by Cwkbj in the expression for Qr,r. 
In this way we obtain 

Qr.r= (1 - ^) Q. ( lllf^lf) (1 - ^corr COs(^/(2 + q)) 
\^corr 



\nq/{2 + q)) 
r(2/(2 + g)) 



where 



T[[l + q)/{2 + q))r 



Ac;/) 



(113) 



= B, 



sin(7r/(2 + q)) 

2a;Nj2/(2+'J) r((3 + g)/(2 + g)) 
r((l + g)/(2 + g))' 



m 



(114) 



see equation ( 1107) , and 
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i^inK/2)- 



'2(2 + ,) /^ 3(4+,) 



Expressions ( 111311151 determine the contribution of the radiative 
region close to the convective base to the overlap integral. Note that 
Acorr <C 1 when u is small. It thus represents a correction to the 
leading term. 

6.4.2 Contribution to the overlap integral from the convective 
envelope 

In order to find this contribution we first determine the constant Cc 
that scales the solution l |98l l, applicable in the convective envelope, 
in terms of the constant C2 that scales the solution in the radiative 
region. This can be done by applying the condition that the radial 
displacement ^ be continuous at r = Vc. Using ( |98t and l llOIt 
evaluated at x = to obtain the displacements just outside and just 
inside r = re respectively, we find that the continuity condition 
gives 



Cc = C2(A'/2a;)-5+i 



r((l + g)/(2 + g))C' 



(116) 



where ( = {r^'^ /pc)xc^^^ F{bi,b2,b3,Xc) and Xc — I — re. 
After using ( 1116) to substitute for Cc in J98t we substitute the re- 
sult into J80b . Making use of ( |92| (. ( 11071 ) and (TTTJ we evaluate the 
integral to obtain 

30, 




^)(1- AeorrCOs(^— 



))Qb 



where 



with 



/c = 



7rp;(A'/2) 2(2 + 5) I^OJ 2(2 + 9) 

r((l + g)/(2 + g))r?/'Ai/4 



J' drr^+'^xl^-' F(b^MM,x) 



s°--2)^vr^ 



CwKB.l; (117) 



(118) 



n 'Xc' ' F(bi,b2,b3,Xc) 

The integral /c(A) is plotted as a function of A in Fig.|5] It will 
be seen that Ic depends rather weakly on A in the considered 
range, being approximately equal to its value for a non-rotating star, 
namely 1.45 • 10"^. 

Note that when A is assumed to be constant (as in the non- 
rotating case) and g = 1, from ( 1117) it follows that Qr (x ojT. 
This scaling was obtained by Zahn (1970) using a different formal- 
ism. 

6.4.3 Analytic expressions for the normalised overlap integrals 

The normalised overlap integrals entering the expressions for the 
tidal energy and angular momentum transfer given by ( 150) are re- 
lated to Qr, a, and n, through 



Qj = aC 



where Qr ~ — ;=• 



(119) 



Note that a is defined through equation ( 170) and its form for the 
dominant mode {i — 1) is shown in Fig.[T] The quantities Q,. can 
be specified as Qr — Qr.r + Qr.c, where Qr,r and Qr,c are given 
in equations ( 1113) and ( 1117) . respectively and we recall that n is 
given by dllO) . 

In order to evaluate these expressions, along with a, we re- 
quire A, both of which are shown in Fig. [T] Be, and Ic, the lat- 
ter two quantities being plotted in Fig. \5\ In addition to these we 




Figure 7. The dependence of Qr on ui for the non-rotating CD model with 
different fitting models for the radial dependence of the Brunt - Vaisala 
presented in table[T] See the text for a description of the different curves. 



require the stellar model dependent quantities, /, A, q, re, and pc. 
These are specified below and in table[T] The above information en- 
ables straightforward calculation of analytic estimates for the nor- 
malised overlap integrals for rotating stars, which match well the 
numerically determined ones, for forcing frequencies less than the 
critical rotation period. 

The values of the integral I are taken to be equal to 14.91 and 
5.898 for the CD and IR models, respectively. The density pc, in 
standard units, at the base of convection zone is equal to 0.1048 in 
the CD model and to 0.1375 in the IR model. It is assumed that the 
base of the convection zone is situated at re — 0.729 in all cases. 



6.5 Quantitative comparison of the the overlap integrals 
obtained numerically and analytically 

The comparison of the analytic estimates for the overlap integrals 
determined from equation ( 1119) with values obtained numerically 
is illustrated in Figures |7llll 

In Figs |7] and [8] we show Qr as a function of cj for the non- 
rotating CD and IR models, respectively. For non-rotating models 
we straightforwardly have Q — Qr- The solid curves give numer- 
ically evaluated overlap integrals with the uppermost one corre- 
sponding to the IR model. For the non-rotating models of the star 
the Cowling approximation is relaxed. The dashed, dotted and dot 
dashed curves give our analytically determined overlap integrals for 
the CD2 and IR2 fitting models (dashed curves), for the CDl and 
IR2 fitting models (dotted curves) and for the CDL and IRL mod- 
els (dot dashed curves), respectively. As seen from this Fig. when 
(J is small our analytic results are an excellent agreement with the 
numerically determined ones, especially for the CD2 and IR2 fit- 
ting models, which will be used later in other Figures. For the CD2 
model the disagreement is smaller than or of the order of 10 per 
cent when oj ^ 0.5 and smaller than 30 per cent when (j ^ 1. For 
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Figure 8. As in Fig.|2]but for the IR model. 



Figure 10. As in Fig. [T] but showing Qr for a rotating star with Qr = 
±0.42. All curves are for the CD model with the CD2 fitting model for the 
radial dependence of the Brant - Vaisala frequency. 
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Figure 9. As in Fig. [T] but showing the contribution of different terms to 
Qr- All curves are for the CD model with the CD2 fitting model for the 
radial dependence of the Brant - Vaisala frequency. 



Figure 11. As in Fig.fTOlbut for the IR model of the star with the IR2 fitting 
model for the radial dependence Brant - Vaisala frequency. 



the IR model the disagreement is even smaller being of the order of 
or smaller than 10 per cent when oj ^ 1.5. 

It is very important to stress that the IR model, describing a 
young zero age main sequence star of one solar mass, has larger 
overlap integrals. This means that tidal interactions are stronger in 
that case. Our analytic approach to the calculation of overlap inte- 
grals allows us to explain this as follows. One can check that the 
main difference in the values of the overlap integrals is determined 



by the the integral / entering in the expression for the norm ( IllOV 
As mentioned above, we have 7 = 14.91 and 5.83 for the CD and 
IR models, respectively. Our analytically determined overlap inte- 
grals are inversely proportional to square root of /. Therefore, the 
overlap integrals of the IR models are larger than the ones of the CD 
models by a factor approximately given by a/14.91/5.83 ~ 1.6. 
Since the energy and angular momentum exchanges through tidal 
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interactions are proportional to square of the overlap integrals, in 
the case of the zero age MS star, these are approximately 2.55 times 
larger in comparison with those obtained for the present day solar 
model. 

In Fig|9]we show different contributions to Qr, for the CD2 
model, separately. The solid and dashed curves are as in Fig.[71 The 
dot dashed curve shows the analytically determined Qr, but with 
with Q,.,c set to zero in ( 1119) . The dotted curve shows Qr, but with 
Qr,r = in ( |119b . Thus, the dot dashed and dotted curves show 
the separate contributions of the radiative region close to the base 
of convection zone and the contribution from convection zone i to 
Qr- As seen from Fig.|9l Qr,r is more important for relatively large 
values of u), while the Qr,c gives a contribution that becomes in- 
creasingly important in the limit of oj —> 0. However, for the range 
of uj shown, both contributions must be taken into account in or- 
der to obtain good agreement between the numerical and analytical 
results. Additionally, we show the overlap integrals of a star mod- 
elled as an n = 3 polytrope, plotted using the dot dot dashed curve. 
Unlike the overlap integrals for the the realistic star, the polytropic 
ones decrease exponentially with decreasing u. Therefore, when g 
modes of sufficiently high order determine the tidal exchange of en- 
ergy and angular momentum, this is much smaller for a polytropic 
star. 

In Fig. [To] we compare numerical and analytic results for a 
rotating CD2 model. In order to delineate the effects of rotation, 
the rotational frequency is taken to be large with \Q,r\ ~ 0.42. 
The solid, dashed, and dot dashed curve show the numerical re- 
sults corresponding to Qr = 0, —0.42 and 0.42, respectively. The 
dotted curves correspond to the analytic calculations and are clos- 
est to their corresponding numerical curves. One can see that the 
agreement between analytical and numerical results is again very 
good for u < 1 for the non-rotating star and the star with prograde 
rotation. The analytic and numerical curves can hardly be distin- 
guished. For the case of retrograde rotation the curves disagree by 
up to 30 per cent in the range shown. The curves corresponding 
to prograde rotation are qualitatively similar to the curves corre- 
sponding to a non-rotating star, the only difference being that they 
indicate somewhat larger values of Qr for a given uj. The curve 
corresponding to retrograde rotation behaves in a non-monotonic 
manner. It has a strong minimum at a; = 0.44, which corresponds 
to 1/ = —1.9. From Fig. [T] it follows that A = 30 for this value 
of v and from equation dlllt it follows that the factor h{r) en- 
tering our analytic expression for Q,. is then equal to zero. This 
explains the non-monotonic character of the behaviour of Qr for 
retrograde rotation. In general, in this case the radial overlap in- 
tegral is smaller than for the non-rotating star. This as well as the 
decrease of the 'angular overlap integral' a for retrograde rotation 
(see Fig.[T]l suppresses the tidal transfer of energy and angular mo- 
mentum for sufficiently large retrograde rotation. Note, however, 
that for retrograde rotation there are other factors, such as a match 
between forcing and mode frequencies that tends to occur for lower 
order modes in comparison to the prograde case, that may amplify 
the transfer, see e.g. Ivanov & Papaloizou (2011) . In Fig. I HI the 
corresponding results are given for the rotating IR2 model. The re- 
sults are qualitatively similar to the previous case. 

In Figs [12] and [T3] we show the dependence of the overlap in- 
tegral defined through ( 1119) (the index j being suppressed) with uj 
set equal to 1 in l lll3t and ( I117t on v, for the CD2 and IR2 mod- 
els, respectively. Note that the overlap integral given by ( |119l > in- 
cludes the factor contributed by the angular overlap integral a{i'). 
The solid curves show Q, the dashed curves show Q calculated 
with the contribution from the convective envelope, Qr,c defined in 
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Figure 12. The value of the overlap integral ni9t for a fixed eigenfre- 
quency lu = 1 plotted as a function of u together with the contributions 
to it from the regions just inside the radiative core and the convective en- 
velope. The CD2 fitting model is used. See the text for a description of the 
different curves. 
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Figure 13. As for Fig.[T2]but for the IR2 model. 



( 11171 ) set to zero and the dot-dashed curves show Q calculated with 
the contribution from the radiative interior, Qr.r defined in ( 11131 ) 
set to zero. One can see from Figsll2landll3lthat at negative values 
of V, Q becomes quite small. The minimum at v ~= 1.9 is de- 
termined by the fact that A ~ 30 there and the integral is strongly 
suppressed as discussed above. 

As follows from equations ( I113t and dllTt the contributions to 
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the overlap integral coming from the radiative interior near base of 
convection zone and the convective envelope are, to a good approx- 
imation, proportional to powers of cj. This can be used to extend the 
results presented in Figs.ll2landll3lto other values of ui by making 
appropriate scaling. Of course, when using such results, one has to 
take into account that the values of lj, the mode order and Q are 
linked together through equation ( 1106) and also the definition of v. 



7 TIDAL EVOLUTION TIMESCALES FOR A BINARY 
WITH A SMALL ECCENTRICITY AND WITH 
ORBITAL AND SPIN ANGULAR MOMENTA 
ALIGNED IN THE MODERATELY LARGE 
DISSIPATION REGIME 

As a simple illustration of application of our formalism we calcu- 
late the timescales for the evolution of the orbital elements of a 
binary system induced by tidal interaction. We consider a binary 
consisting of a primary star of mass M, on which tides are exerted 
by the secondary of mass A/2 that is approximated as being a point 
mass. For simplicity, we assume below that the orbit has a small 
eccentricity , e, and the stellar rotation frequency is small in com- 
parison with the characteristic frequency fJ* . 

We shall assume that the action of dissipation on the eigen- 
modes is sufficiently strong for the system to be in the regime of 
moderately strong viscosity discussed above. Then we can use ( |51| l 
to calculate the orbital energy dissipation rate due to tides. In the 
case of Sun-like stars considered here, this may require the excited 
waves to be of large enough amplitude to break as they propagate 
towards the stellar centre (see Barker 2011 and references therein 
for a discussion). Although such a condition may not be satisfied, 
assuming that it is should provide the most optimistic estimate of 
the potential long term tidal evolution due to the excitation of these 
modes. Note too that the action of turbulence in convection zones is 
likely to be ineffective for the short orbital frequencies we consider 
on account of its long characteristic time scale. 

This problem has been discussed in a similar setting by e.g. 
Zahn (1977) and Goodman & Dickson (1998). However, they did 
not include the effects of Coriolis forces, taking account of rota- 
tion only through a shift in tidal forcing frequency. Our formalism 
includes the effects of Coriolis forces under the traditional approx- 
imation and so allows us to assess their significance. Moreover, the 
formalism of Zahn (1977) has been mainly applied to more massive 
stars having a convective core and radiative envelope in contrast to 
the solar type stars explicitly considered here (see also Papaloizou 
& Savonije 1997). Although our general approach is applicable to 
this case, the asymptotic calculation of the overlap integrals needs 
to be considered separately. 

In sections |2]- [4] the stellar angular velocity was taken to be 
positive so that prograde normal modes had positive frequencies 
while retrograde normal modes had negativefrequencies. However, 
for the discussion of noiTnal modes and overlap integrals from sec- 
tion|4]onwards, it was convenient to exploit the symmetry that en- 
abled retrograde modes to have a positive frequency, with the sign 
of il being reversed to become negative. In this section we consider 
the tidal interactions of a star with a definite positive angular veloc- 
ity, accordingly from now on we revert to the convention of section 
|2]-|4]for which retrograde modes have negative frequencies. 



7.1 Equations for the evolution of the semi-major axis and 
eccentricity 

Equations governing the evolution of the semi-major axis a, and 
the eccentricity, e, follow from consideration of the conservation of 
energy and angular momentum of the system. We begin by noting 
that tidal interaction with M* induces rates of change of energy, 
Ei, and angular momentum, Lc, for the orbit given by equation 
dSlb . Conservation of energy and angular momentum implies that 
these produce changes in a and e according to 



GfiM'ia 



^Ei, 



GfiM^ee 



= -E/(l 77) ^ 0,orbLc 



(120) 



2a2 " a ''* 2 

Here the eccentricity is assumed to be small so that powers of e 
higher than second order may be neglected and fi = M2/M, is the 
mass ratio. In what follows, we find it convenient to rewrite these 
equations in the form 

d 2 e 1 



T ' 



where 



Te^ 



aEe 



T, 



Ta^ 



GfiM'i 
aEi 



and 



Ee = Elil — ) — Q.orbLa 



(121) 



(122) 



(123) 



When the rotation frequency ^ is specified we can use equa- 
tions (JSTJ and the expressions for the coefficients Am,k given in 
Appendix A to estimate the time scales Ta and Te. We consider be- 
low two important cases 1) the case of synchronous rotation Q, — 
Qorb and 2) the case of a non-rotating primary 57 = 0, labelling 
the corresponding time scales by indices (1) and (2), respectively. 
Since we shall use dimensionless units, we write Q = \/Mt.Rt,Q, 
where Q is the dimensionless overlap integral, we recall that we 
have suppressed the index j. 

7.2 The case of synchronous rotation 

Here we have fl — Qorb- As indicated in Appendix A it turns out 
that to leading order in e, terms of higher order in e than linear in 
the expansion of the tidal potential may be neglected when substi- 
tuting into the series ^. Thus, the term proportional e^/2 in l ll23t 
results in a higher order contribution and can be neglected. We thus 
have Ee ~ Ec. The (m, k) pairs that remain to be considered are 
(0, 1), (2, 1) and (2,3). The conditions for resonance associated 
with a normal mode ( I41t associated with these are all included in 
the expression cj = ujm,k ~ k^lorb — rn^ = zt^lorb- The corre- 
sponding values of the parameter u = ±2. As seen from Figs.[T] 
and |2] when either v = —2 and m — 2 01 v = ±2 and m = 
we have A ~ 34. This leads to a strong suppression of the corre- 
sponding overlap integrals determined by the factor h{r) oc 1 — ^, 
see the discussion above, equation illlt and the discussion of Fig. 
llOlfor the m — 2 case. Therefore, their respective contributions to 
the energy and angular momentum exchange rates can be neglected 
and we take into account below only the dominant (2, 3) term. 

Equations l lSlt and ( 1123b then give to the leading order in ec- 
centricity 



E, = -- 



/TttV /GMaV 
"40 V aO 



M,Ri 



\duj/dj\ 



and 



El = 3Ee 



(124) 



where functions of eigenfrequencies uj enclosed in square brackets 
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Figure 14. The evolution time scale T^ ' calculated for the CD and IR 
models compared to the result of Goodman & Dickson (1998), for /i = 1. 
See the text for the description of different curves. 



Figure 15. The tidal Q' parameter for orbital circularisation. The solid 
and dashed curves are respectively for CD2 and IR2 models in the case of 
synchronous rotation, the dotted and dot-dashed curves are respectively for 
the CD2 and IR2 models with no rotation. 



[..]fc are assumed to be evaluated for cj — kQort and m — 2 from 
now on. 

The term [|da;/dj|]fc can be evaluated setting n = j in equa- 
tion ( 1106b and adopting the limit n 3> 1, noting that A varies only 
weakly with z^ in the region of interest being close to 4.5 ( see 
Fig.[T]), with the result that 



\dLu/dj\]k = 






(125) 



where the integral I is expressed in natural units. Substituting ( 1125) 
in ( I124t we obtain 

Ei = ^^ [Q'] ^ , and E. = —— [q'] ^ , 



40 



where 



E, = VAI—^, T 



"orb" 



(126) 



Now we use ( I126t to obtain the tidal evolution time scales given by 
( II22I 1 thus obtaining 



n(l) 



n{l) 



40T, 



1477r 



[Q-1.- 



= -e~"T, 



2^(1) 



where 



(l+/i)'/' /Porb^V^'o-l 



(127) 



(128) 



VA^.1 

and the orbital period Porb = 27rfi^,,'^. Thus, for this configuration 
we have Ta 2> Te when e ^ 1. 

In order to calculate the overlap integral entering (1127) for 
m = 2 we use the results of preceding section setting values of 
parameters determining the overlap integral to correspond to ;/ = 



2. Explicitly, we use ai ~ 1, A = 4.5, le 
and Ic = 1.4 • 10"^ 

We show the evolution timescales T, 



(1) 



-0.1155, Be 



calculated for the CD 



and IR models in Fig.[T4] The short dashed, dotted, long dashed and 
dot dashed curves are for the CD2, CDL, IR2 and IRL fitting mod- 
els. One can see from Fig.ll4lthat there is little difference between 
the curves corresponding to different fitting models. The solid line 
plots a minor modification of the expression of Goodman & Dick- 
son (1998) who give T^° = 8 ■ 10^ Pj^;,, with Port in days, for a 
system of two tidally interacting Sun-like stars of equal mass. 

Since we assume that the tides are raised only on one compo- 
nent, we plot Tf ^ = 1.6 ■ 10*PX-(,//i in Fig.[T4]and others below. 
Our results are seen to give larger values of Te for a given Port , e.g. 
for Porb ~ 4 days, the results differ by a factor of five to ten. Thus 
we have Te « 2.7 • 10** yr, 1.15 ■ lO'^yr and 2.35 • 10^ yr from the 
modified Goodman & Dickson (1998) expression and for our cal- 
culations corresponding to the CD and IR models, respectively. We 
remark that Ogilvie & Lin (2004) indicate a possible reason why 
the Goodman & Dickson expression gives such an underestimate 
for Te and we recall that the effects of Coriolis forces are not in- 
cluded. Note that the IR model gives larger values of Te irrespective 
of the fact that its overlap integrals are larger than those of the CD 
model, for a given forcing frequency. The is due to a smaller radius 
of the early Sun compared to its present value. The younger star is 
more compact and, accordingly, less susceptible to tidal influence 
from a perturber. 

Finally let us recall that for large values of Port, the appropri- 
ate normalised overlap integrals corresponding to the CDL and IRL 
fitting models are proportional to P^^^ . From equation ( I127t it 
follows that in this case Te oc Pj^-i, as for the Goodman & Dickson 
(1998) expression. 

Finally we calculate the effective tidal Q' for the star by com- 
paring our circularisation rate to that of Goldreich & Soter (1966). 
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Figure 16. The evolution time scales T„ ' and T^ ' calculated for the CD 

and IR models compared the modified expression of Goodman &; Dickson 

(1998). See the text for the description of the different curves. 



(2) 
Figure 17. The evolution time scales T,, plotted for q 



10- 



10"'' and q 
lively. 



10 using the solid, dashed and dotted curves, respec- 



This comparison leads to the expression 

285 






(129) 



987r%/A/[Q2]i(l + ^)i/2 

We plot Q as a function of orbital period for \i = \m Fig. [T5] 
This quantity is often assumed to be a constant. However, Fig. 1151 
shows that this varies between lO" and 10* as the orbital period 
varies from 1 — 10 days. 

7.3 The case of non-rotating primary 

Calculation of the time scales Ti and Ti, proceeds analogously 
to the previous case. Since fJ = 0, the forcing frequencies are of 

the form LOm.k ~ kQ.orb- 

f21 

In order to calculate T^ we can set e = in the expressions 
for the potential expansion coefficients given in Appendix A. Then 
only the A2.2 term with corresponding forcing frequency 012,2 = 
2Q,orb contributes to the series. The calculation is straightforward 
with the result that 



Ei 



SnE:, 



[Q' 



40 L-- J2' 

and, accordingly. 



n(2) 



4or. 



[Q- 



(130) 



(131) 



It is convenient to separate the calculation of Ee, and, accord- 
(2) 
ingly, Te into two stages. First we evaluate Ei — QortLc- This 

quantity is determined by terms in the potential expansion with 
forcing frequencies uJo,i = ^orb, W2,i = ^orb and lj2,3 = SQ,orb- 
The respective coefficients Am,k are given in Appendix A. Sub- 
stituting these into ( 15 It , noting that in the non rotating case, the 
overlap integrals are independent of m, we obtain 

Ej \lorbi-'c ^ 
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Figure 18. Same as Fig.[T7]but T^ ' is shown. 
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20 



Q' 



|da;/dj| 



497r^ 



0" 



\du)/dj\ 



(132) 



Note that the term associated with the forcing frequency Sflort 
gives the leading contribution since the normalised overlap inte- 
grals strongly decrease as a; — >■ 0. 

Now we calculate e^_E//2. To leading order in e^ it is suffi- 
cient to take into account only the term in the potential expansion 
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with forcing frequency 0^2,2 = 2Q.orh, thus from equation ( I130t we 
obtain 

Stt^ /GM2' 



y-- 



0' 



da;/dj| 
From ( I123t and J132t we then find that 



(133) 



£e 



where 



-M''"'-" 



2/GM2 

V a3 



) F{^or 



(134) 
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\ ^" 1 


-3 
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[\d^/dj\\ 


2 


r Q' 1 

[\du/dj\\ 


[id^MlJ 



(135) 



Now we make use of il25t . applicable in this case, to ehminate the 
frequency derivatives with respect to j in l ll35b finally obtaining 



where 



J'i^orb) — 



and, accordingly. 



{St^^: 



-AQ' 



+Ym.}' 



n{2) 



20 



n> = — T,/j^(n„,6). 



(136) 
(137) 

(138) 

In Fig. [16] we show the characteristic times Ti and T^ , 
for fJ. — 1. The short dashed and dotted curves give Te for 

the CD2 and IR2 models respectively. The long dashed and dot 

(2) 
dashed curves give Ta , respectively for the CD2 and IR2 models. 

The solid curves correspond to the modified Goodman & Dickson 

(2) (2) 

(1998) expression. We see that Ta is much larger than T^ ' m all 
cases. This is obviously due to the higher forcing frequency SClorb 
effective in the latter case. It is interesting to note that the numeri- 

(2) 

cal values of Te are either smaller than, or close to the result of 
Goodman & Dickson (1998) regardless of the fact that our results 
were calculated for a non-rotating primary, while the Goodman & 
Dickson (1998) result corresponds to a primary in a state of syn- 
chronous rotation. 

In Figs.ll7landll8lwe show the evolution timescales T^ and 
Ta for mass ratios appropriate for exoplanetary systems. The solid, 
dashed and dotted curves correspond to M2 = Mj, O.lMj and 
lOMj. As seen from these plots, dynamical tides in the regime 
of moderately large viscosity are potentially efficient only when 
orbital periods are rather small. Thus, the eccentricity is damped in 
a time less than lO^yrs only for Port < 1.85, 2.56 and 3.48daj/s, 
for M2 = O.lAfj, Mj and lOMj, respectively. Similarly the semi- 
major axis can decay in a time less than l(fyrs only for Port < 
1.24, 1.69 and 2.3days respectively. 

We obtain Q' for the star by comparing our circularisation 
rates e/e to that of Goldreich & Soter (1966). In this case we find 



Q 
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/ 27r 
[Portn 



987rVA/[02]i(l + ^)l/2 \Porb^, 

for the case of synchronous rotation and 

855 



Q = 



47rjr(no,6)(l + ^)i/2 






(139) 



(140) 



in the non-rotating case. We plot Q as a function of orbital period 
for /i = 1 in Fig. [TS] This quantity is found to vary between 10* 
and 10^ as the orbital period varies from 1 — 10 days. 

Finally we make a rough estimate of the effect of inertial 
waves that could be excited in the convective envelope in the case of 



synchronous rotation. We note that Ogilvie & Lin (2007) calculate 
the Q' associated with inertial waves, for a solar model with a ro- 
tation period of ten days as a function of tidal forcing frequency as 
viewed in the rotating frame, appropriate to forcing with a m = 2 
quadrupole potential (see their Fig.3). For the relevant forcing fre- 
quencies id = ±Q,orb, Q' ^ W , implying that the effect of iner- 
tial modes may be comparable to those due to rotationally modified 
gravity modes at the longest orbital periods considered ~ 10 days. 
But note that if the tidal dissipation due to these waves scales as fi 
as indicated in Ogilvie (2013), their effect may be somewhat less 
important at shorter orbital periods. 



8 CONCLUSIONS AND DISCUSSION 

In this paper we have calculated the general tidal response of a uni- 
foiTnly rotating star and used it to find the rate of transfer of energy 
and angular momentum between the star and the orbit. The results 
are expressed as a sum of contributions from normal modes each 
being proportional to the square of the associated overlap integral. 

We considered the response due to an identifiable regular 
dense spectrum of normal modes such as one provided by the low 
frequency rotationally modified gravity modes. We obtained ex- 
pressions ( I5U for the energy and angular momentum transfer rates 
appropriate to the moderately viscous case for which the waves as- 
sociated with the modes are damped before reaching an appropriate 
boundary (the centre for a radiative core and the surface for a radia- 
tive envelope) or understood in the averaged sense as described in 
section U.5.11 In that case they are independent of the details of the 
dissipation process. In addition we obtained expressions l |52b appli- 
cable to the case of very weak dissipation where resonant responses 
may occur. We also gave corresponding expressions ( 154) that can 
be used to evaluate the energy and angular momentum exchanged 
as a result of a flyby in a parabolic orbit. These can be applied to ei- 
ther the tidal capture problem or the problem of tidal circularisation 
at large eccentricity (see Ivanov & Papaloizou 2011). 

Our expressions may be applied to either Sun-like stars with 
radiative cores and convective envelopes or more massive stars with 
convective cores and radiative envelopes. However, in this paper 
we considered only Sun-like stars in detail. We developed expres- 
sions for the overlap integrals in the low frequency limit, under the 
traditional approximation. These were evaluated analytically fol- 
lowing a WKBJ approach that yielded eigenfrequencies and eigen- 
functions for the rotationally modified gravity modes. Good agree- 
ment was found between analytic estimates of the overlap integrals, 
obtained from equation ( |119t making use of equations ( 1113) and 
( 1117b . and values obtained through direct numerical integration. 
These showed that for prograde forcing the overlap integrals in- 
creased as the rotation rate of the star increased and accordingly 
tidal interaction strengthened. The opposite, but more pronounced 
trend, was found for retrograde forcing (Figs 110113] . Knowledge 
of the overlap integrals and mode spectrum enables calculation of 
tidal energy and angular momentum exchange rates with the orbit 
as a function of tidal forcing frequency in the moderately viscous 
regime. 

For Sun-like stars with inner radiative regions, our formal- 
ism allows one to find the tidal response of low frequency gravity 
modes without numerical evaluation of eigenfrequencies and over- 
lap integrals, and thus find the orbital evolution of tidally interact- 
ing stars of this type. In order to implement it, it suffices to know 
the behaviour of the Brunt - Vaisala frequency close to the base of 
convection zone, the radius and the value of the density at that loca- 
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tion, as well as the integrals I ~ J,," ^N, and, for rotating stars, 
the quantities determining the angular properties of the overlap in- 
tegrals and given in Fig[T] In general, it is also necessary to know 
the decay rate of the WKBJ due to either dissipation or non-linear 
processes. Provided this information is given as a function of stel- 
lar age, one can construct more realistic models of tidal interaction, 
that allow for finite values of the eccentricity, a non zero inclination 
of the stellar rotation axis to the orbital plane , and stellar evolution 
by semi-analytic methods. It is important to generalise the analytic 
evaluation of the overlap integrals to apply to stars with different 
structure, eg. more massive stars with convective cores. This will 
be undertaken in future work. 

It is interesting to note that in our formalism, the dependence 
of energy and angular momentum transfer rates on semi-major axis 
is determined by the dependence of the normalised overlap in- 
tegrals Q on the forcing frequency co. The standard dependence 
(Zahn 1977) is recovered only when Q oc aj^^'*", which may, in 
general, not be satisfied. Although in the models we consider in 
this paper, the dependence of Q on oj is very close to a power law 
with index 17/6 in the limit oj — >■ 0, it differs from this for short 
period orbits. It may also differ for stellar models with a different 
radial structure. This emphasises the necessity of analytical and nu- 
merical studies of the behaviour of the overlap integrals for stellar 
models with different masses and ages. 

We used our formalism to determine the tidal evolution 
timescales for a binary consisting of a primary Sun-like star in 
which the tidal response was excited in the moderately large dis- 
sipation regime. Assuming synchronous rotation and that inwardly 
propagating rotationally modified gravity waves are dissipated 
close to the centre, the quantity, tJ^' ((l + /x)/2)^/^/^, where Ti^' 
is the circularisation time, was found to vary between 10'^ — 10^'^ yr. 
as the orbital period varied through the range 1 — 5d. The corre- 
sponding value of Q' defined by Goldreich & Soter (1966) varied 
through the range 3 x 10^ - 10*. 

Similarly when the primary was not rotating, TJ ((1 + 
n)/2)^'^ /fj,, Tg being the circularisation time in this case, was 
found to vary between 10* — lO^'yr. as the orbital period varied 
through the range 1 — 5d. The quantity Q' has a similar range 
of variation as is found in the synchronous case. In this non syn- 
chronous case, the quantity Ta ((1 + fi)/2)^'^ /fi where the time 
scale for inspiral is Ta varied between 10® — 10^" yr as the orbital 
period varied through the range 1 — 5d. Thus the latter process is 
unlikely to have the possibility of being effective for masses in the 
planetary regime (fi < lO^'^) and orbital periods exceeding 2.5d. 

In this paper when tackling the problem with periodic forcing 
potential we assume that the mode energy doesn't grow with time 
and the amount of energy supplied by tidal interaction is equal to 
the amount of dissipated energy. This condition may however be 
broken either very close to the resonances or in the case when the 
distance between two eigen frequencies is smaller than the inverse 
of a characteristic orbital evolution time. In the latter case the ap- 
proach based on the Fourier transform may be more adequate as 
discussed in eg. Ivanov & Papaloizou (2004a). 
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APPENDIX A: THE COEFFICffiNTS A^.k FOR A 
COPLANAR SLIGHTLY ECCENTRIC ORBIT 

Equation ^ gives the perturbing tidal potential acting on the pri- 
mary component in the form 



* = r^ ^ (A„,fee-'"'".^'Fr + cc) 



(Al) 



An expansion of the gravitational potential due to the perturbing 
secondary in spherical harmonics yields the coefficients Am,k- We 
consider the case when the eccentricity of the binary orbit is small, 
such that powers of e higher than the second may be neglected, 
and when the rotational axis is perpendicular to the orbital plane. A 
convenient expression for the tidal potential that can be used to de- 
termine the Am,k under these conditions can for example be found 
in Zahn (1977). The non-zero coefficients are found to be 
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(A2) 



Here M2 is the mass of the secondary and a is the orbital semi- 
major axis. Note that for the purpose of evaluating the complex 
conjugate in ^ we have A^^ j. = Am,k = A-rn,-k- We also re- 
mark that in order to determine the rate of change of the semi- major 
axis correct to zero order in eccentricity and the rate of change of 
eccentricity correct to first order in eccentricity, terms oc e^ may be 
dropped when working with these coefficients. 



APPENDIX B: COMPARISON OF ASYMPTOTIC 
EXPRESSIONS FOR TORQUES AND ENERGY 
DISSIPATION RATES GIVEN BY NORMAL MODE SUMS 
WITH THOSE OBTAINED USING RADIATION 
BOUNDARY CONDITIONS 

Here we consider the response of a spherically symmetric stel- 
lar model with rotation under the traditional approximation. The 
governing equation for the linear response to a forcing potential 

*m,fe exp(— iaj„,ftt + irrKJ)) is 

-^i,fc^m,fe-2ia;m,fc(n-f)fx^„^fc+C^^_fc = -V*„,fc,(Bl) 



This follows from equation (|5} provided $7 is replaced by its radial 
component $7 f . Making this replacement amounts to making the 
traditional approximation described in section |4] Of course when 
n is set to zero, the problem reduces to the problem of finding the 
response of a spherically symmetric non rotating star. Thus con- 
sidering equation ( IB It enables us to consider both rotating stars 
with the adoption of the traditional approximation and non rotating 
spherical stars. 

The action of the operators C and V respectively on ^^ j. and 
^m,k is as indicated by equation l |15l l. It is taken as read that in 
addition to time dependence through a factor Gxp(— iojm.fct), each 
of ^j7i,fe ey.p{ira(f>) and ^^ f. exp(im0) have separable angular de- 
pendencies related to the Hough functions for the specified m and 
with an appropriate value of A (see section|4j. These dependencies 
will not be indicated below. We go on to discuss the response due 
to a single forcing term. 

When the response has been obtained for a range of assumed 
angular dependencies corresponding to related values of A, a total 
response can be found using the orthogonality of the Hough func- 
tions (see section |4!4t . Total torques or energy dissipation rates are 
then obtained by summing contributions associated with each value 
of A involved in the representation of the perturbing potential in 
terms of Hough functions. Note also that the values of A that are 
expected to dominate are associated with the potentially resonant 
mode spectrum, namely that associated with rotationally modified 
gravity modes, or simply g modes in the case of a non rotating star. 

For simplicity of notation, the subscripts ra,k on all quanti- 
ties will be suppressed from now on. In addition we adopt spheri- 
cal polar coordinates (r, 9, cf)) for the purpose of discussion in this 
appendix. 



Bl Torque evaluated with radiation boundary conditions 

The equation to be solved is equation JBU . Noting the angular de- 
pendence indicated above, A solution can be found by the method 
of variation of parameters where we write 



^e, ^ = ^(r)^i + B(r)^2 



(B2) 



In the above ^^ is a solution of the homogeneous form of equation 
JBU {^ set to be zero) that satisfies the correct boundary conditions 
enforcing regularity at r = and ^2 is a solution of the homoge- 
neous form of equation dBlt that satisfies the appropriate boundary 
conditions at the surface r ~> Rs- Here we take the latter solution 
to have to correspond to an outgoing wave. In this case physical 
dissipation is supposed to take place for r ^^ Rs- But note that we 
expect a corresponding parallel analysis to hold when ^-^ is taken 
to correspond to an ingoing wave while regularity conditions are 
applied to ^2 for '' ^^ Rs- In that case physical dissipation is sup- 
posed to take place for r ~ 0. In any case the function A{r) is cho- 
sen to vanish for r — >■ _Rs and B{r) is chosen to vanish for r — )■ 0. 
In addition the functions A{r) and B{r) satisfy the constraint 



^i-V4(r) + ^2-VB(r) = 0. 



(B3) 



This ensures that derivatives of A and B do not appear when the 
pressure perturbation P' is evaluated (see below). 

We remark that for the case when dissipation occurs near the 
surface considered here, when uj corresponds to a free oscillation 
eigenvalue, ^^ will correspond to the free normal mode. As we 
consider the asymptotic low frequency limit, where the spectrum 
is dense, and we assume a degree smoothness in the quantities of 
interest, ^-^ may be assumed to correspond to the normal mode 
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with eigenfrequency closest to the forcing frequency. Thus for com- 
parison with the modal summation approach described above, we 
identify ^j with ^ . We further comment that the equilibrium tide 
is implicitly included in the decomposition l lB2b . But the impor- 
tant contribution for computing the net torque is the outward going 
wave component and the equilibrium tide does not affect this. 

We make use of the expressions for C^, P'{$,) and A'^^ given 
by equations ( 155b . l l56b and ( I57l > and we take the adiabatic exponent 
Fi to be constant, for convenience, as in the main text. It follows 
from these that for any pair of solutions of the homogeneous form 
of equation JBlb such as ^^ and $2 the latter equation implies that 



14^12 = / {p'{i2)il - p'{ai2) ■ ds, 



(B4) 



where Sr is an interior spherical surface of radius r, is constant. In 
addition for the solution of the inhomogeneous equation of interest, 
we have 

/ {P'imi - P'iCiH) dS- f pil- V*dT = 0.(B5) 

The expression ( IB4b is associated with the conservation of wave 
action for freely propagating disturbances. 

From equation JBlb we find that the total torque communi- 
cated to the forced body is given by 



To(^) = -mXm 



P'{^)e ■ dS 



(B6) 



Using equation l lB6b with the fact that the forced solution l lB2t be- 
comes 5^2 ^s the surface of the star is approached, the angular 
momentum transfer rate to oscillatory disturbances as a result of 
the tidal interaction may be evaluated as 

To(^) = \B{R,)\'To{i^) (B7) 

Furthermore, using I IB51 we can express B{Rs) in the form 



B(Rs 



W12 



Thus the tidally induced torque may be written as 



To(0 = 



|Wl2|2 



(B8) 



(B9) 



To relate this to expressions given in the main text, we now use the 
WKBJ approximation to find alternative expressions for some of 
the quantities in ( 



WKBJ solutions 

As follows from equation ( |87b in the main text after neglecting the 
1 as compared to A^ /a; in the last bracket on the left hand side, in 
the small uj^ limit, free oscillations are governed by the following 
equation for ii — r^P^'^^S^r, Cr being the radial component of the 
displacement, 

d / n dv'' 
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= -X « 



p2/ri rfj. \p2/ri ^j. 
where 
2 _ TV'Ap" 

^ ^2^2p4/ri • 

We remark that the pressure disturbance is given by 
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(B12) 



APi/ri dr 

As follows from equation ( 188b . away from turning points, the 
WKBJ solution for vi — r^P^'^^ £,r,i, where ^r,i is the radial 
component of ^ j , can be written in the form 



Vl 



W2 



(B13) 



where $ — N\^/{rij), rt is the inner boundary radius of the 
wave propagation zone and 771 is a constant phase determined by 
the boundary condition at r = 0. 

The WKBJ solution for V2 = r^P^^^^ ^r,2, where ^^,2 is the 
radial component of the outgoing wave solution ^2 j takes the form 



expj 



V2 



(-x: 



$dr + r)2 



X 



1/2 



(B14) 



where 1)2 is another constant phase factor. We may use the above 
WKBJ solutions to evaluate the constant quantities To ($2) ^ti^ 
W\2 in the WKBJ approximation in the form 



To{i2) = 



and 



(B15) 



■[¥12 = —r—cxp 1(772 — r;i — 7r/2). (B16) 

We remark that here we have assumed without loss of generality 
that the integration of the square of the angular factors over the 
spherical polar angles is normalised to unity. Using the above ex- 
pressions in ( |B9b we obtain 



To(^) 



u]\lo\ 



(B17) 



The WKBJ spectrum 

To proceed further we recall that from the free and undamped form 
ofdBB 

u^Uj - C^j - t^j 6^i=0 (B18) 

We use the fact that in the WKBJ limit the spectrum is given in the 
large j limit by the phase integral relation 



\/A/^"Afr-idr 



JTT 



= ^o/j, 



(B19) 



where the ± alternative corresponds to the positive and negative 
frequency modes respectively, and see equation ( |103b of the main 
text with n = j > 0. Thus, noting that A is frequency dependent, 
we have 



, / ,. 1 , ^3 9A 



2 / 



(B20) 



In order to re-express the above expression in terms of the mode 
norm, we begin by noting that if we replace A'^ by A'^(l + e) in 
( |B19b and regard loj to be a function of e, we have 



d^./de|,=„(l-g-^)=a;,|.=o. 



(B2I) 



However, the same result as ( IB21t can be obtained from dBlSb by 
multiplying by p^* and integrating over the fluid volume and then 
treating the introduction of e from the point of view of perturbation 
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theory. Noting that the operators are self adjoint, we find from this 
that 



2a;, 



pN'^li, ■ r\^dT 



(B22) 



Using the WKBJ solution <B13I ) to evaluate the right hand side of 
l lB22b . identifying ui and ojj with u)j„ and ^j^ with ^ and making 
use of l |B20l (. | |B2T] | and l lBl2l l we obtain 



du)j/dj\ja 



27V,„A 



(B23) 



where we recall the definitions of the inner product (|4] and norm 
(|26}. 



Using ( |B23l l together with l |B17t enables the latter to be re- 
duced to the form 

|2 



To(0 = ^ 






(B24) 



Recalling again the definition of the inner product this can be writ- 
ten in the alternative form 



^ 2^^mc.,„($;jVM/)^ 
The rate of energy dissipation is then 



Ec = -u,^To{i)/m = 



iVjo|dcjj(j)/djjj=j, 



(B25) 



(B26) 



This should be multiplied by a factor of two to account for the 
response to the complex conjugate of the applied forcing term as 
specified by for example by equation (|6}. This is seen from the form 
of l |B6l ) to give an identical contribution. The result then becomes 
identical to a term corresponding to single values of m and k in the 
summation given by ( 145b . 
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